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•  This  report  considers  the  fire  environment  that  would  result  from  a 
megaton-yield  nuclear  weapon  explosion.  An  analysis  (developed  in  Part  L)  that 
treats  the  physics  of  the  burning  zone  and  the  volume  immediately  above  it 
(turning  region)  is  used  to  predict  the  velocity,  temperature,  and  pressure 
fields  of  large  area  fires. 

A  sensitivity  study  explores  the  influence  of  turbulence,  radiation,  fire 
size,  and  burning  intensity  on  the  mean  temperature  levels  and  velocity  fields. 
The  results  show  hurricane-force  velocities  developing  as  the  fire  size  or 
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burning  rate  is  increased.  A  sample  calculation  illustrates  the  change  in  fire- 
wind  velocities  as  the  fire  evolves  over  time. 

Calculations  of  the  burning  region  for  three  model  urban  areas  show  the 
influence  of  building  density  and  urban  sprawl  on  the  resulting  fire  environment. 
An  additional  set  of  predictions  accounts  for  reduction  of  the  fire  intensity  by 
blast  in  the  urban  center.  For  the  latter  cases,  the  temperature  distribution 
is  changed  markedly,  though  the  magnitude  of  the  induced  fire  winds  is  not 
appreciably  reduced. 
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PREFACE 


This  report  is  the  second  of  the  two-part  documentation  of 
Pacific-Sierra  Research  Corporation's  analysis  of  the  large  urban 
fire  environment.  This  part  presents  calculations  of  the  character¬ 
istics  of  nuclear-weapon-ignited  fires.  Part  I  develops  the  theory 
underlying  the  analysis.  All  work  was  performed  for  the  Federal 
Emergency  Management  Agency  under  contract  EMW-C-0747;  the  technical 
monitor  was  Dr.  David  Bensen. 

The  authors  gratefully  acknowledge  the  contributions  of  Dr. 
Harold  L.  Brode. 
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SUMMARY 


This  report  considers  the  large-fire  environment  that  would  occur 
in  an  urban  area  subject  to  a  nuclear  weapon  explosion.  The  effects 
of  system  parameters  are  explored  in  a  sensitivity  study,  and  results 
for  three  model  cities  are  presented. 

The  examples  are  characterized  by  extensive  areas  simultaneously 
burning,  strong  buoyancy,  and  large  temperature  gradients.  Several 
such  fires  occurred  during  World  War  II.  Though  those  fires  were 
dramatic  in  intensity  and  destructiveness,  each  involved  a  relatively 
small  area.  A  nuclear  weapon  explosion  could  generate  a  far  larger 
area  fire  and  a  more  severe  fire  environment.  This  report  is  intended 
to  define  such  large  fires.  The  results  should  be  applicable  for 
damage  evaluation,  formulation  of  shelter  requirements,  and  rescue 
planning. 

The  calculations  are  based  on  the  theory  developed  in  Part  I  of 
this  report,  which  is  applicable  to  the  fire  zone  and  the  volume 
immediately  above  it  (turning  region).  The  effects  of  variable  area 
heating,  turbulence,  strong  buoyancy,  large  temperature  changes,  and 
radiation  are  treated.  The  induced  fire  winds  and  rapid  temperature 
changes  at  the  fire  periphery  are  uniquely  determined  by  the  use  of 
jump  conditions.  Simulations  of  the  Hamburg  firestorm  and  a  large 
Flambeau  fire  agreed  well  with  available  data. 

The  parametric  analysis  considers  a  large  area  fire  and  the  ef¬ 
fects  of  fire  sii,^  ueating  rates,  mixing  coefficients,  and  hot  gas/ 
smoke  radiation.  The  results  show  the  influence  of  those  variables 
on  the  induced  fire  winds,  mean  temperature,  and  pressure  gradients. 

In  general,  an  increase  in  either  the  fire  size  or  heating  rate  raises 
the  mean  temperature  levels  and  the  induced  fire-wind  velocities. 

For  the  larger  heat  release  rates  or  fire  sizes,  the  attendant  in¬ 
creases  in  mean  temperature  and  velocity  are  limited  by  compressi¬ 
bility  effects. 

Fires  such  as  may  result  from  a  megaton-yield  explosion  are 
analyzed  for  three  model  urban  areas.  Each  city  is  characterized  by 
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a  high-density  center,  a  surrounding  belt  of  mixed  residential/ 
industrial  construction,  and  a  lower  density  suburban  belt.  Each 
model  city  portrays  a  different  degree  of  building  density  and  urban 
3prawl.  The  results  illustrate  how  a  particular  city  geometry  affects 
the  velocity  and  temperature  fields  for  a  given  fire.  An  additional 
series  of  computations  considers  reduction  of  the  fire  area  by  severe 
blast  damage  and  debris  formation.  For  those  calculations,  complete 
burning  was  allowed  in  an  annular  area,  with  the  fire  intensity  signif¬ 
icantly  reduced  in  the  center. 

Finally,  the  model  is  employed  to  estimate  the  behavior  of  the 
velocity  and  temperature  fields  as  a  function  of  fire  evolution. 

Those  calculations  may  indicate  the  most  appropriate  periods  for 
effecting  rescue  operations  as  well  as  provide  an  estimate  of  time- 
dependent  shelter  loadings. 
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Cp  *  specific  heat  capacity  at  constant  pressure 

Eg  =  average  heat  released  per  unit  weight  of  all  combustibles 

=  effective  kinematic  viscosity  (for  turbulent  flow) 

fg  =  building  density  (ratio  of  area  covered  by  buildings  to 
total  area) 

F  =  dimensionless  measure  of  QH  [defined  by  Eq.  (13)] 

g  =  gravitational  acceleration 

H  =  maximum  height  of  flames 
* 

k  =  reciprocal  of  graybody  radiation  mean  free  path 
kp  =  effective  thermal  conductivity  (for  turbulent  flow) 

Kp  =  dimensionless  heat-diffusion  coefficient 
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Qa  *  areal  heat-addition  rate 
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r  =  radial  position  coordinate  (dimensionless,  except  in 
Figs.  24  through  26) 

R  »  fire  radius 

T  »  temperature  (dimensionless) 

=■  ground-level  atmospheric  temperature  in  far  field 

T  *  maximum  cemperature 
max 

u  =  radial  velocity  (dimensionless) 

u  =  maximum  radial  velocity 
max 

U  =  radial  velocity  scale 

v  =  vertical  velocity  (dimensionless) 

v  =  maximum  vertical  velocity 
max 

y  =  vertical  position  coordinate  (dimensionless) 

Y  “  specific  heat  ratio 

2 

5  *  dimensionless  constant,  U  /(P  /p  ) 

cl  3 

tp  =  radial  dependence  of  model-city  heat  release  rates  in 
blast-modified  cases 

p  *  density  (dimensionless) 

Pa  *  ground-level  atmospheric  density  in  far  field 

a  =  radiation  coefficient  (dimensionless  constant), 
4TraT'!tk*/Q 

cl 

a  =  Stefan’s  constant 

symbols'^  used  in  documentation  of  prediction  algorithm 

*  intermediate  variable  [defined  by  Eq.  (A. 37)] 

F,  G,  H  *  vectors  of  functions  that  are  to  be  driven  to  zero  (by 
proper  choice  of  x)  in  Newton  iterations,  and  that 
represent  discretized  model  equations  on  individual 
lines  of  constant  y 

T>lathematical  variables  only — not  FORTRAN  labels. 
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L. _ INTRODUCTION 


This  report  presents  predictions  of  the  temperatures,  pressures, 
and  high-speed  winds  created  by  large  urban  fires .  The  dependence  of 
those  quantities  on  fire  size,  burning  rate,  and  various  other  param¬ 
eters  is  explored,  and  fires  in  model  U.S.  cities  are  examined. 
Simulations  in  which  nuclear-weapon-ignited  fires  are  extinguished 
in  the  center  by  blast  are  compared  with  those  in  which  the  fires  con¬ 
tinue  to  burn.  The  model  developed  in  Part  I  of  this  report  was  used 
to  calculate  all  values.  That  model  may  be  readily  extended  to  obtain 
estimates  of  oxygen  depletion  ard  noxious  gas  buildup  as  well  and, 
hence,  a  fairly  complete  baseline  description  of  the  environment  facing 
civil  defense  personnel  when  large  urban  areas  are  burning. 

Most  previous  research  concerning  the  hydro thermodynamics  of 
free-burning  fires  has  been  restricted  to  the  weakly  buoyant  free- 
convection  plume.  The  resulting  theories  describe  the  basic  flow  in 
the  middle  and  upper  parts  of  a  long,  thin  plume  over  a  small  fire, 
and  may  be  relevant  to  some  portion  of  the  convection  column  generated 
over  a  large  urban  fire  [e.g.,  Morton,  Taylor,  and  Turner,  1956;  Murgai 
and  Emmons,  1960;  Yokol,  1960;  Nielsen  and  Tao,  1965].  However,  as  dis¬ 
cussed  in  Part  I,  such  theories  are  inapplicable  in  and  around  the  com¬ 
bustion  zone,  so  they  cannot  be  used  to  predict  the  environment  and 
high-speed  surface  winds  induced  by  large  fires. 

Smith,  Morton,  and  Leslie  [1975]  present  global  computations  of 
the  hydrothermodynamic  environment  near  the  combustion  zone  as  well 
as  above  it.  In  their  calculations,  all  components  of  the  fire¬ 
generated  flow  field — e.g.,  surface  inflow,  the  convection  column, 
the  far  field — are  considered  collectively,  and  the  strong  coupling 
between  dynamic  forces  and  the  induced  surface  winds  is  demonstrated. 

In  and  around  the  combustion  zone,  however,  density  changes  are 
assumed  to  be  small  (Boussinesq  approximation),  heat  is  input  at  the 
boundary  rather  than  in  volume,  and  not  much  resolution  is  sought. 


Further,  the  computations  are  not  specifically  designed  for  the 
analysis  of  large  urban  fires. 

In  a  more  recent  study,  Brode,  Larson,  and  Small  [1982]  adopt 
the  basic  computational  approach  of  Smith,  Morton,  and  Leslie  in 
studying  the  mesoscale  motions  generated  by  large  fires,  but  eliminate 
some  of  the  stated  deficiencies.  Their  work  is  directly  aimed  at  the 
definition  of  the  large  urban  fire  environment,  and  arbitrary  changes 
in  temperature  and  density  are  allowed  in  the  near-fire  region.  The 
resolution  in  that  region  is  still  rather  coarse,  however. 

The  component  analysis  introduced  by  Small,  Larson,  and  Brode 
[1981]  provides  an  alternative  to  the  global-calculation  analysis  in 
which  fine  resolution  near  the  fire  can  be  obtained  economically. 

Using  that  approach,  the  components  of  the  flow  field  are  analyzed 
separately,  and  then  combined  by  means  of  appropriate  matching  condi¬ 
tions.  The  special  features  of  each  component  are  thus  considered  in 
greater  detail,  and  resolution  is  enhanced. 

The  key  element  in  the  component  analysis  is  the  description  of 
the  environment  in  and  around  the  fire,  since  the  entire  flow  is 
driven  by  the  hydrothermodynamic  interactions  in  that  region.  The 
turning-region  model  developed  in  Part  I  provides  the  key  description. 
That  model  decouples  from  those  for  the  other  flow  components,  and 
thus  can  be  used  independently  to  make  detailed  predictions  of  the 
near-fire  environment.  Such  a  capability  extends  our  understanding 
of  fire  dynamics  and  is  the  basis  for  the  present  analysis. 

We  use  the  turning-region  model  to  generate  the  large  urban  fire 
predictions.  In  that  model,  a  volume  heat  source  represents  the  net 
effect  of  the  combustion  process,  and  the  induced  flow  is  taken  to  be 
axisymmetric  and  quasi-steady.  The  flow  is  also  assumed  compressible, 
to  permit  arbitrary  changes  in  temperature  and  density.  A  one- 
parameter  eddy-viscosity  model  is  used  to  describe  the  turbulent 
transfer  of  heat  and  momentum,  and  a  graybody  approximation  to  model 
hot  gas  and  smoke  radiation.  Finally,  jump  conditions  are  derived 
to  account  for  the  rapid  changes  in  physical  quantities  at  the  fire 
periphery.  They  effect  model  problem  closure,  and  decouple  the  turning- 
region  analysis  from  that  of  other  component  flows.  An  iterative 
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finite-difference  scheme,  documented  in  the  Appendix,  is  employed  to 
solve  the  resultant  boundary  value  problem. 

We  study  the  large  urban  fire  environment  in  two  ways — by  varying 
parameters  and  by  simulating  three  model  cities.  In  the  parametric 
analysis,  excursions  are  made  about  a  baseline  in  which  the  fire  radius 
is  taken  to  be  10  km  and  the  characteristic  heat  release  rate  is  con¬ 
sistent  with  that  of  the  1943  Hamburg  firestorm  [DCPA,  1973].  Predic¬ 
tions  of  hurricane-force  winds  are  typical  in  the  parameter  excursions 
as  well  as  the  model  city  simulations,  which  consider  fires  of  radius 
12  km  with  spatially  varying  heating  rates  on  the  order  of  those  in  the 
parametric  analysis.  That  is  believed  to  be  about  the  size  of  the  fire 
that  would  result  from  a  near-surface  nuclear  burst  of  about  1  Mt 
[Johnson  and  Larson,  1982]. 

Section  II  presents  the  results  of  the  parametric  analysis.  As 
expected,  temperatures  and  fire-wind  velocities  increase  with  an  in¬ 
crease  in  either  fire  size  or  burning  rate,  though  the  fire-wind 
variations  are  not  as  rapid  as  the  linear  scaling  law  of  Part  I  sug¬ 
gests.  Also  as  expected,  temperatures  and  winds  increase  with  a 
decrease  in  radiation  intensity;  winds  also  increase  with  a  decrease 
in  the  magnitude  of  the  eddy  coefficient  for  turbulent  momentum  trans¬ 
fer.  The  temperature  field  is  relatively  insensitive  to  changes  in 
that  coefficient,  however,  and  both  winds  and  temperature  are  rela¬ 
tively  insensitive  to  changes  in  the  corresponding  coefficient  for 
turbulent  heat  transfer. 

The  basic  flow  pattern  also  seems  to  be  relatively  insensitive 
to  variations  in  the  specific  dependence  of  the  heat  release  rate  on 
spatial  position,  though  such  variations  modify  velocities  slightly 
and  may  significantly  change  the  temperature  field.  That  behavior 
is  found  when  model  predictions  for  a  uniformly  heated,  circular  fire 
are  compared  with  those  for  an  annular  fire.  The  comparison  indi¬ 
cates  how  partial  blast  extinguishment  would  modify  the  hydrothermo¬ 
dynamic  environment  generated  by  a  nuclear-weapon-ignited  fire.  It 
also  provides  a  first  look  at  the  catastrophic  effect  of  multiple 
weapon  bursts. 


« 


Section  III  describes  the  model  city  simulations.  Three  model 
cities  are  defined:  one  that  is  lightly  built-up,  intended  to  repre¬ 
sent  new,  sprawling  cities;  one  that  is  heavily  built-up,  intended  to 
represent  old,  congested  cities;  and  one  of  intermediate  building 
density.  For  each  city,  spatially  dependent  fuel  loadings  and  peak- 
period  burning  rates  for  a  baseline  fire  are  estimated,  and  the  result¬ 
ing  temperatures,  pressures,  and  fire  winds  predicted.  Those  calcu¬ 
lations  are  repeated  for  a  scenario  in  which  the  fire  is  ignited  by  a 
nuclear  burst  and  the  city  partially  destroyed  by  the  attendant  blast. 

Somewhat  surprisingly,  despite  a  significant  difference  in  model- 
city  temperature  predictions  for  the  "blast"  and  "no  blast"  simula¬ 
tions,  the  corresponding  difference  in  induced  fire  winds  is  fairly 
small.  Such  behavior  suggests  that  the  winds  and  wind  damage  result¬ 
ing  from  nuclear-weapon-ignited  fires  may  be  relatively  independent 
of  the  degree  of  blast  extinguishment — unless,  of  course,  extinguish¬ 
ment  is  nearly  complete. 

As  expected,  however,  temperatures  and  fire-wind  velocities  all 
increase  with  building  density.  For  a  given  fire  size,  the  large-fire 
threat  facing  civil  defense  workers  will  thus  be  more  severe  in  the 
tall,  heavily  loaded  cities  than  in  the  short,  lightly  loaded  cities. 

In  general,  however,  the  shorter  cities  sprawl  out  over  greater  land 
areas  than  do  caller  ones  of  comparable  population,  and  are  thus  capable 
of  supporting  larger  fires.  The  large-fire  threat  for  any  given  popu¬ 
lation  size  may  thus  be  greatest  in  some  of  the  shorter  cities,  espe¬ 
cially  if  the  fires  result  from  multiple  nuclear  bursts.  Comparisons 
between  specific  cities  must  be  made  individually. 

Finally,  Sec.  IV  presents  a  sample  application  of  the  quasi-steady 
model  of  Part  I  to  obtain  time-dependent  predictions.  The  sample  case 
considers  multiple-fuel-bed  Flambeau  fires,  the  largest  area  fires  for 
which  a  reasonable  body  of  technical  data  exists  [Countryman,  1969; 
Palmer,  1981 J . 
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II.  PARAMETRIC  ANALYSIS  OF  LARGE-FIRE  ENVIRONMENT 


Using  Che  curning-region  model  of  Part  I,  this  section  explores 
the  dependence  of  the  large-fire  environment  on  fire  size  (height  and 
radius),  burning  rate  (spatial  average  and  form),  degree  of  turbulent 
mixing,  and  degree  of  hot-gas  radiation. 

Briefly,  with  r,  y,  u,  v,  p,  T,  and  P  denoting  dimensionless 
radial  and  vertical  position,  radial  and  vertical  velocity,  density, 
temperature,  and  perturbation  pressure,  respectively,  the  turning- 
region  equations  are  as  follows: 


~  (rpu)  +  ^  (rpv)  =  0  , 


(la) 


/  9u  9u\  3P  /  1  9  /  9u  \  u  \ 

plu  37  +  v  3?/- - 3?  +  MlU¥F(r  3F/“7j  * 


+  Ap  =  0  , 


p  (u  If +  v!f) =  q(r’ y)  - a(T*  - 1}  +  Kl(f  £(r  If))  ’ 


(lb) 


(lc) 


(Id) 


PT  -  1  ,  (le) 

where  the  dimensionless  parameters  are 
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and 


M, 


V°» 

UR 


kl/cPpa 

UR 


(3) 


Here,  R  and  H  denote  fire  radius  and  maximum  flame  height;  Q  the 
scale  and  q(r,  y)  the  spatial  distribution  for  the  rate  of  volumetric 
heat  addition;  and  P^,  pa,  and  Ta  ground-level  atmospheric  pressure, 
density,  and  temperature  in  the  far  field,  respectively.  Further, 

yv  * 

a  is  Stefan  s  constant,  k  is  the  reciprocal  of  the  graybody  radiation 

mean  free  path  (assumed  constant),  and  u>  and  k  are  dimensional  eddy 

1  *  + 

coefficients  for  turbulent  momentum  and  heat  transfer. 

We  consider  the  dependence  of  model  predictions  on  the  size  param¬ 
eters  R  and  H,  the  burning  rate  parameters  QH^*  and  q(r,  y) ,  the  radia¬ 
tion  mean  free  path  k*  ,  and  the  dimensionless  eddy  coefficients 

and  K  .  The  terms  £,/p  and  k./c^p  remain  constant,  except  when 
X  x  a  X  a  a 

variations  with  and  are  specifically  studied. 

BASELINE  ANALYSIS 

As  a  baseline  case,  we  consider  an  axisymmetric  fire  with  the 
following  parameters: 


R  =■ 

10 

km  , 

(4a) 

H  - 

100 

m  , 

(4b) 

QH  =* 

5.7 

/  •y 

*  10^  cal/m  -sec  , 

(4c) 

As  in  Part  I,  g  is  the  acceleration  due  to  gravity,  y  the  spe¬ 
cific  heat  ratio,  and  U  the  radial  velocity  scale.  The  vertical 
velocity  is  scaled  with  (H/R)U,  temperature  and  density  with  Ta  and 
p_,  and  perturbation  pressure  with  6P  ,  where  0  «  U^/(P  /p  ). 


a  ‘  ‘  a'  a  '  a' 

’QH  is  used  as  a  parameter  instead  of  Q  because  areal  heat  re¬ 
lease  rates  are  encountered  much  more  frequently  than  volumetric  rates 
in  the  large-fire  literature. 


4 


—  /  — 


(1.6 

for 

0  £  y  £  0.25 

q (r ,  y)  =  ■.  1.6 

(f  (1  -  y)) 

for 

0.25  sy  s.1.0  , 

(  0 

for 

y  a  1.0 

Mx  =  Kx  -  0.2  ,  (4e) 

*-l 

k  =20  m  .  (4f ) 

The  heat  release  race  in  this  case  is  the  same  as  chat  in  our  simula¬ 
tion  in  Part  I  of  the  1943  Hamburg  firestorm,  and  is  characteristic 
of  the  rates  estimated  in  Sec.  Ill  for  a  variety  of  U.S.  cities.  The 
parameter  QH  is  the  nominal  specific  scale  rate  suggested  by  DCPA  [1973]. 
The  form  of  q(r,  y)  represents  maximum  heating  in  and  around  the  fuel 
zone  and  a  decrease  in  heating  with  increased  altitude.  All  other 
parameter  choices  used  in  Eqs.  (4)  are  of  the  same  order  as  those  used 
in  the  model  city  analysis. 

From  Eqs.  (2)  and  (4),  the  radial  velocity  scale  for  the  baseline 
case  is 


U  =  67.2  m/sec  , 


(5) 


and 


A  =  0.217  , 

C  =  0.110  .  (6) 

The  temperatures,  pressures,  and  velocities  predicted  for  this  case 
are  summarized  in  the  vector,  contour,  and  profile  plots  of  Figs.  1 
through  5.  The  flow  field  is  illustrated  in  Fig.  1.  Temperature  and 
pressure  contours  are  plotted  in  Figs.  2  and  3,  respectively.  The 
temperature  attains  a  maximum  in  the  fire  center  at  the  top  of  the 
maximum  heating  zone  (y  =  0.25).  It  then  decreases  rapidly  with 
increased  altitude,  the  flow  becoming  weakly  buoyant  above  three 


r 


Altitude  (m  x  100) 
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Horizontai  valocrty  Im/wcl  :  -33.6  0.0 


Radius  (km  x  10) 

Fig.  4--Radial  velocity  profiles  for  baseline  fire 


Radius  (km  x  10) 


Mg.  o— Vertical  velocity  profiles  for  baseline  fire 
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flame  heights.  As  expected,  the  pressure  drop  is  maximum  at  ground 
level.  Profiles  of  the  horizontal  inflow  induced  by  the  pressure  drop 
are  presented  in  Fig.  4,  and  those  of  the  resulting  upflow  in  Fig.  5. 
The  maximum  inflow  occurs  on  the  ground  at  the  fire  periphery.  The 
upflow  profiles  are  nearly  "top-hat." 

In  most  parameter  excursions,  the  hydrothermodynamic  predictions 
are  qualitatively  the  same  as  those  for  the  baseline  case.  The  re¬ 
sults  are  summarized  below. 

DEPENDENCE  ON  FIRE  SIZE  AND  BURNING  RATE  SCALE 

We  consider  first  the  dependence  of  the  large-fire  environment 
on  R,  H,  and  QH — the  parameters  that  collectively  describe  the  fire. 

For  nuclear-weapon-ignited  fires,  the  effective  value  of  R  will  depend 
on  the  number  and  yields  of  bursts,  ambient  atmospheric  conditions, 
various  characteristics  of  the  preblast  urban  area  (e.g.,  flammability 
of  exposed  combustibles,  susceptibility  of  structures  to  secondary 
fires),  and  a  variety  of  lesser  factors  [Johnson  and  Larson,  1982]). 
Parameters  QH  and  H  will  depend  on  the  type  and  distribution  of  com¬ 
bustibles  in  the  postblast  environment,  as  detailed  in  Sec.  III. 

Tables  1  through  3  summarize  the  variations  that  we  consider  in 
fire  size  and  burning  rate,  and  the  corresponding  changes  in  model 
parameters  that  follow  from  Eqs.  (2)  and  (3).  In  most  excursions, 
the  variations  are  of  the  type  in  Table  2,  where  the  parameter  under 
study  is  simply  doubled  or  halved.  The  additional  variations  in 
Tables  1  and  3  are  considered  in  order  to  refine  the  excursions  that 
we  believe  to  be  of  greatest  intrinsic  interest — i.e.,  those  defining 
the  dependence  of  fire  winds  and  temperatures  on  fire  radius  (weapon 
coverage)  and  intensity  (concentration  of  combustibles). 

Figures  6  through  11  summarize  the  results  of  the  fire  size  and 
burning  rate  excursions.  The  dependence  of  fire  winds  and  tempera¬ 
tures  on  fire  size  and  intensity  is  plotted  in  Figs.  6  through  8.  As 
expected,  the  maximum  induced  velocity  umax  and  the  maximum  temper¬ 
ature  T  both  increase  with  either  radius  or  intensity.  The  in¬ 
max 

creases  are  nearly  linear  for  relatively  small  radii  and  heating  rates, 
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Table  1 

PARAMETER  DEPENDENCE  ON  FIRE  RADIUS 


R  (km) 

U  (m/sec) 

A 

a 

V  K1 

2.3 

16.8 

3.4720 

0.110 

3.200 

5.0 

33.6 

0.8680 

0.110 

0.800 

10.0 

67.2 

0.2170 

0.110 

0.200 

15.0 

100.8 

0.0964 

0.110 

0.089 

20.0 

134.4 

0.0534 

0.110 

0.050 

Table  2 

PARAMETER  DEPENDENCE  ON  FIRE  HEIGHT 


H  (m) 

U  (m/sec) 

A 

a 

V  K1 

50 

134.4 

0.0272 

0.055 

0.4 

100 

67.2 

0.2170 

0.110 

0.2 

200 

33.6 

1.7360 

0.220 

0.1 

Table  3 

PARAMETER  DEPENDENCE  ON  BURNING  RATE  SCALE 


QH-, 

(kcal/m— sec) 

U  (m/sec) 

A 

a 

V  K1 

14.35 

16.8 

3.4720 

0.4400 

0.80 

28.70 

33.6 

0.8680 

0.2200 

0.40 

57.40 

67.2 

0.2170 

0.1100 

0.20 

86.10 

100.8 

0.0964 

0.0825 

0.15 

114.80 

134.4 

0.0534 

0.0550 

0.10 

Fire  radius  (km) 


Fig.  6— Dependence  of  maximum  radial  velocity  and  temperature 

on  fire  radius 


"'"max 


Burning  rate  scale  (kcal/m^  sec) 


Fig.  7--Dependence  of  maximum  radial  velocity  and  temperature 

on  burning  rate  scale 


Fire  height  (m) 

Fig.  8— Dependence  of  maximum  radial  velocity  and  temperature 

on  fire  height 


but  tail  off  markedly  for  fires  of  much  greater  size  and  intensity. 

In  fact,  u  and  T  appear  to  asymptote  toward  finite  values  as 

R  or  QH  -*■  00 .  Surprisingly,  the  apparent  asymptotic  values  of  umay 

as  either  R  or  QH  +  “  are  almost  identical.  T  ,  however,  exhibits 

max 

a  stronger  dependence  on  QH  than  on  R. 

Figure  8  describes  the  variation  of  u  and  T  with  fire  height 

max  max 

As  H  increases,  u  increases  but  T  decreases.  Such  behavior  is 
max  max 

presumably  explained  as  follows.  With  QH  fixed,  Q  decreases  when  H 
increases.  Temperature  decreases  correspondingly,  because  of  the 
shorter  characteristic  residence  times  of  fluid  particles  in  the  heat¬ 
ing  zone.  Lower  temperatures  might  be  expected  to  decrease  velocities 
as  well.  However,  a  smaller  fraction  of  the  constant  total  heat  re¬ 
lease  QH  is  radiated  away  at  lower  temperatures,  so  higher  kinetic 
energies  can  be  supported. 

The  results  in  Fig.  8  must  be  interpreted  with  some  care.  When¬ 
ever  H  is  varied,  QH  may  be  varied  as  well.  For  example,  when  build¬ 
ings  have  the  same  construction  and  contents,  both  H  and  QH  should 


Maximum  perturbation  pressure  (psi)  Maximum  perturbation  pressure  (psi) 


Maximum  perturbation  pressure  (psi) 


200  H 


200  H 


Fire  height  (m) 

Fig.  11— Dependence  of  maximum  perturbation  pressure  and  vertical 

velocity  on  fire  height 


increase  with  the  number  of  stories.  In  such  cases,  the  results  in 
Fig.  7  may  be  more  applicable. 

Figures  9  through  11  plot  the  dependence  of  maximum  pressure 

drop  P  and  vertical  velocity  v  on  R,  QH,  and  H.  As  expected, 
max  J  max  r 

vertical  velocities  generally  increase  with  radial  velocities. 
Greater  pressure  gradients  foster  greater  radial  velocities,  and  are 
less  strongly  correlated  with  temperature  variations  than  Eq.  (lc) 
suggests. 


DEPENDENCE  ON  BURNING  RATE  SPATIAL  DISTRIBUTION 

We  now  briefly  explore  the  effects  of  different  spatial  depen¬ 
dences  of  the  heat  release  rate  q(r,  y) .  We  consider  several  sinu¬ 
soidal  and  annular  variations  in  radial  dependence,  but  just  one  in 
vertical  dependence.  Further  variations  are  considered  in  the  model 
city  analysis  in  Sec.  III. 

Table  4  summarizes  the  variations  in  q(r,  y)  and  the  resulting 
nydrothermodynamic  changes.  As  expected,  a  relatively  high-frequency 
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Table  4 

VELOCITY,  TEMPERATURE,  AND  PERTURBATION  PRESSURE  MAXIMA 
FOR  VARIOUS  BURNING  RATES 


(q(r,  y)  \  u  TP  v 

t~, —  ~t  {  ^ - J  max  max  max  max 

"  r*  ^  'baseline/  (m/sec)  (°K)  (psi)  (m/s ec) 


Baseline 

1.0 

37.2 

5  70 

0.156 

3.26 

Sinusoidal 

Low-freq uency 

1.0 

-  COS 

47rr 

38.1 

588 

0.160 

3.50 

High-frequency 

1.0 

-  cos 

20iTr 

37.2 

573 

0.157 

3.28 

Annular 

o  o 

for  0 
for  0. 

£  r  < 
5  s  r 

0.5) 

*  1) 

33.9 

515 

0.119 

2.46 

radial  variation  on  the  baseline  heat-release  distribution  has  little 
effect.  Surprisingly,  however,  a  fairly  low-frequency  radial  varia¬ 
tion  also  produces  trivial  changes  in  velocity,  temperature,  and  per¬ 
turbation  pressure.  Those  results  suggest  that  the  large-fire 
environment  may  be  more  sensitive  to  the  basic  scale  for  heat  release 
than  to  the  specific  spatial  distribution  of  the  heat  release  rate 
(see  Fig.  7) . 

That  conclusion  is  further  supported  by  a  comparison  of  the  re¬ 
sults  in  Table  4  for  the  annular  fire  with  those  for  the  circular 
baseline  fire.  The  ground  area  covered  by  the  annular  fire — and  hence 
its  total  heat  release — is  25  percent  less  than  that  of  the  baseline 
fire.  The  corresponding  drop  in  maximum  temperature  difference  from 
ambient  (i.e.,  in  T  -  T  )  is  just  slightly  less — 20  percent.  The 
drop  in  ^max  is  significantly  less — only  9  percent.  Such  behavior 
suggests  that  the  maximum  winds  and  temperatures  generated  by  nuclear- 
weapon-ignited  urban  fires  may  be  relatively  insensitive  to  changes 
in  the  geometry  and  loading  of  the  blast-damaged  region,  provided  of 
course  that  the  blast  does  not  nearly  extinguish  the  fire. 

Figures  12  and  13  compare  the  temperature  and  flow  fields  of  the 
baseline  and  annular  fires.  Since  the  fire  is  located  only  in  the 
annulus,  the  highest  temperatures  occur  near  the  inner  edge  of  the  fire 
annulus  (r  =  0.5),  rather  than  on  the  symmetry  axis.  The  resulting 
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Fig.  13  Flow  fields  of  baseline  and  annular  fires 


temperature  field  is  thus  qualitatively  different  from  that  of  the 
baseline  fire.  Further,  as  Fig.  13  shows,  winds  in  the  annular  fire 
do  not  blow  toward  the  fire  from  both  edges  of  the  annulus.  Rather, 
they  blow  inward  (toward  the  symmetry  axis)  and  upward  at  all  points, 
as  do  those  in  the  baseline  fire.  By  means  of  heat  convection  and 
branding,  annular  fires  might  therefore  be  expected  to  spread  inward. 

Similarly,  a  cluster  of  separated  large  fires,  such  as  could 
result  from  multiple  nuclear  bursts  over  a  large  city,  might  coalesce 
and  engulf  much  of  the  intervening  region.  The  patterns  of  airflow 
and  fire  spread  envisioned  for  such  situations  are  sketched  in  Figs. 
14  and  15.  The  cluster  of  fires  represented  in  Fig.  14  nearly  forms 
an  annulus,  as  indicated  by  the  dashed  circles.  The  fire-wind  flow 


Fig.  15— Radial  airflow  and  fire  spread  patterns  suggested  for 
cluster  of  three  large  area  fires 


should  thus  be  nearly  axisymmetric  and  of  the  type  illustrated  in 
Fig.  13.  Correspondingly,  the  fire  should  spread  inward,  and  may 
ultimately  cover  much  of  the  area  inside  the  inner  dashed  curve.  As 
shown  in  Fig.  15,  the  flow  and  spread  patterns  for  asymmetric  clusters 
should  be  less  axisymmetric,  but  still  basically  inward.  Such  hy¬ 
potheses  could  presumably  be  verified  by  extensive,  three-dimensional 
numerical  computations. 


We  now  consider  a  variation  of  heat  release  rate  with  altitude. 
Specifically,  we  compare  the  results  for  the  baseline  [q(r,  y)  defined 
by  Eq.  (4)]  with  those  for  the  more  uniform  case  where  q  is  independent 
of  y— and  hence  constant.  That  comparison,  summarized  in  Table  5,  re¬ 
flects  the  results  illustrated  in  Figs.  8  and  11.  For  less  concen¬ 
trated  heat  releases,  temperatures  are  lower,  but  pressure  drops  and 
velocities  are  greater. 


Table  5 

VELOCITY,  TEMPERATURE,  AND  PERTURBATION  PRESSURE  MAXIMA 
FOR  BASELINE  AND  UNIFORM  FIRES 


Type  of 
Fire 

q( 

r,  y) 

u 

max 

(m/sec) 

T 

max 

(°K) 

P 

max 

(psi) 

V 

max 

(m/sec) 

( 

1.6 

for 

0  s  y  s  0.25) 

Baseline  . 

1.6(|  (1  -  y)) 

for 

0.25  sy  s:  1/ 

37.2 

570 

0.156 

3.26 

( 

0 

for 

y  ss  1  * 

Uniform 

1.0 

50.9 

497 

0.231 

7.19 

DEPENDENCE 

ON  TURBULENCE 

AND 

RADIATION 

Finally,  we  consider  the  variation  of  the  large-fire  environment 

with  the  eddy  coefficients  M  ,  K.  and  the  graybody  radiation  param- 

eter  k  .  The  results  are  summarized  in  Figs.  16  through  23. 

The  dependence  of  fire  winds  and  temperatures  on  eddy  coefficients 

is  demonstrated  by  Figs.  16,  17,  and  18.  As  expected,  the  maximum 

radial  velocity  decreases  with  an  increase  in  the  eddy  coefficient  for 

momentum  transfer  (Fig.  16).  Correspondingly,  the  maximum  temperature 

increases.  An  increase  in  the  eddy  coefficient  for  heat  transfer  might 

also  decrease  temperatures,  and  hence  increase  velocities.  In  Fig.  17, 

such  variations  are  shown  to  occur,  but  so  weakly  that  both  T  and 

max 

u  must  be  considered  almost  independent  of  K, .  When  M,  and  K,  are 
max  111 

varied  together,  as  in  Fig.  18,  the  results  are  thus  similar  to  those 
with  only  M^  varied  (cf.  Fig.  16). 


Eddy  viscosity 


Fig.  16— Dependence  of  maximum  radial  velocity  and  temperature 
on  eddy  coefficient  for  momentum  transfer 


f  max 


Heat  transfer  coefficient 


Fig.  17— Dependence  of  maximum  radial  velocity  and  temperature 
on  eddy  coefficient  for  heat  transfer 
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imum  perturbation  pressure 


Eddy  viscosity 


Maximum  perturbation  pressure 


Eddy  coefficients 
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Figure  19  indicates  the  dependence  of  fire  winds  and  temperatures 

* 

on  the  radiation  parameter  k  .  Both  winds  and  temperatures  increase 

with  k*  — as  expected,  since  the  effective  scale  for  radiative  energy 

* 

loss  a  is  proportional  to  k  [see  Eqs.  (1)  and  (2)]. 

The  dependence  of  pressure  drops  and  vertical  velocities  on  M.. , 

*  1 

K^,  and  k  is  shown  by  Figs.  20  through  23.  As  in  previous  comparisons 
(Figs.  6  through  11  and  16  through  19),  variations  in  those  parameters 
are  strongly  correlated  with  variations  in  radial  velocity. 
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III.  MODEL  CITY  ANALYSIS 


Using  the  turning-region  analysis  of  Part  I,  we  obtain  predic¬ 
tions  of  the  large-fire  environment  for  three  model  U.S.  cities. 

City  W  is  lightly  built-up,  and  intended  to  represent  new,  sprawling 
cities.  City  E  is  heavily  built-up,  and  intended  to  represent  old, 
congested  cities.  City  M  is  of  intermediate  building  density.  For 
each  city,  two  cases  are  considered:  a  baseline  fire  and  one  modi¬ 
fied  by  blast.  In  all  cases,  however,  the  fire  radius  is  taken  to 
be  12  km — believed  to  be  about  the  size  of  the  fire  that  would  result 
from  a  1  Mt  near-surface  nuclear  burst  [Johnson  and  Larson,  1982], 

DEFINITION 

Few  metropolitan  areas  are  axisymmetr'ic .  Nevertheless,  most 
cities  have  a  main  business  district  with  high-rise  office  and  apart¬ 
ment  buildings,  surrounded  by  lower  density  tracts.  Our  model  cities 
reflect  that  situation.  Each  has  three  regions:  a  tall  central  city; 
a  residential/industrial  belt  of  intermediate  height  around  the  cen¬ 
tral  city;  and  a  low,  primarily  residential  outer  belt.  We  refer  to 
those  regions  as  the  central  city,  inner  belt,  and  outer  belt, 
respectively . 

Table  6  summarizes  the  geometry  of  the  model  cities,  including 
both  the  horizontal  extent  of  each  city  region,  the  corresponding 
vertical  extent  N^  (average  number  of  stories),  and  the  building 
density  f  .  We  assume  that  the  top  and  bottom  stories  are  each  4  m 

DU  ^ 

high  and  all  other  stories  are  3  m  high.  The  values  for  building 
density  f  in  each  model  city  region  are  comparable  to  those  in  the 
DCPA  manual,  and  thus  considered  representative  of  U.S.  cities  [DCPA, 
1973].  The  values  are  all  less  than  0.45,  which  approximates  the 
densely  packed  Hamburg  districts  devastated  in  the  firestorm  of  1943 
(see  Part  I) . 

^The  average  heights  of  the  outer-belt  buildings  in  cities  W 
and  E  are  taken  to  be  6  and  10  m,  respectively. 
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Table  6 

SIZE  AND  DENSITY  OF  MODEL  CITY  REGIONS 


Region 

City 

W 

M 

E 

Radial 

Dimension  (Tan) 

Central  city 

0-1 

0-1 

0-2 

Inner  belt 

1-3 

1-4 

2-6 

Outer  belt3 

3-12 

4-12 

6-12 

Average  Number  of  Stories,  N 

Central  city 

8 

12 

16 

Inner  belt 

3 

4 

6 

Outer  belt 

1.5 

2 

2.5 

Building  Density, 

Central  city 

0.30 

0.35 

0.40 

Inner  belt 

0.25 

0.30 

0.35 

Outer  belt 

0.15 

0.20 

0.25 

Radius  of 

outer  belt  may  extend 

beyond  12  km. 

but  the  fire  does 

not. 

Calculated 

as  the  ratio  of  land 

area  covered  by  buildings  to  total 
land  area. 


The  corresponding  geometry  of  the  baseline  model-city  fire  simu¬ 
lations  is  summarized  in  Figs.  24  through  26.  In  each  figure,  the 
shaded  area  represents  “.he  fuel  zone  (defined  by  Table  6)  and  the 
hatched  area  represents  the  remainder  of  the  combustion  zone.  The 
total  height  of  the  combustion  zone  in  the  outer  belt  is  assumed  to 
be  five  times  the  fuel  zone  height."^"  Since  turbulence  will  break  the 
flame  envelope  and  tend  to  keep  flame  heights  uniform  [Thomas,  1963], 
the  ratio  of  combustion  zone  height  to  fuel  zone  height  in  the  central 
city  is  taken  as  only  half  as  much;  the  combustion  zone  height  in  the 
inner  belt  is  assumed  to  be  the  average  of  the  other  two. 

x 

'As  noted  in  Part  I,  a  combustion-zone-height  to  fuel-zone-height 
ratio  of  about  five  was  observed  for  the  largest  Flambeau  fires  and 
the  1943  Hamburg  firestorm. 
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The  baseline  areal  heating  rates  (see  Figs.  24  through  26)  for 
the  central  city,  inner  belt,  and  outer  belt  are  estimated  as  follows. 
The  average  fuel  loading  per  story  Lg  in  each  city  region  is  defined 
in  Table  7.  The  corresponding  average  total  areal  fuel  loading  L^  is 
then  computed  from 


L 


T 


L 

s 


N  x 
s 


Bu 


(7) 


with  values  for  N  and  f_  taken  from  Table  6.  The  resulting  total- 
s  Bu 

loading  estimates  are  also  listed  in  Table  7.  The  central  cities  and 
outer  belts  are  somewhat  taller  and  shorter,  respectively,  than  were 
those  in  the  devastated  Hamburg  districts,  and  have  higher  and  lower 
total  loadings.  The  inner  belts  are  defined  to  be  on  the  order  of  the 


Table  7 

FUEL  LOADING  IN  MODEL  CITY  REGIONS 


City 

Region 

w 

M 

E 

Average  Fuel  Loading  per  Story,  Lsa 
( lb/ft2-story ) 

Central  city 

20 

20 

20 

Inner  belt 

18 

18 

18 

Outer  belt 

16 

16 

16 

Total  Areal  Fuel  Loading, 
(lb/ft2) 

lifJJ 

Central  city 

48.0 

84.0 

128.0 

Inner  belt 

13.5 

21.6 

37.8 

Outer  belt 

3.6 

6.4 

10.0 

aBased  on  DCPA  [1973]  estimates, 
which  take  into  account  both  the 
structures  and  their  contents. 


Hamburg  districts,  but  less  built-up.  Similarly,  the  inner-belt  load- 

2 

ings  tend  to  be  somewhat  less  than  the  32  lb/f t  estimated  for  Hamburg 
[DCPA,  1973]. 

Finally,  in  each  model  city  region,  the  fire  is  assumed  to  con¬ 
sume  90  percent  of  the  available  fuel  in  3  hr.  The  areal  heating 
rate  Q,  is  thus  computed  from 

q  »  L  *  /— 90%_\  (g) 

gA  T  \  3  hr  /  B  ’ 

with  L_  taken  from  Table  7  and  E_,  the  average  heat  released  per  unit 

1  D 

weight  of  combustibles,  taken  as  [DCPA,  1973] 

E„  *  8000  Btu/lb  .  (9) 

For  the  blast-modified  fire  simulations  represented  in  Figs.  24 
through  26,  the  weapon  burst  is  assumed  to  occur  over  the  city  center, 
leveling  many  buildings  in  the  central  city  and  inner  belt.  A  cen¬ 
tered,  low-level  burst  of  a  1  Mt  weapon  ould  produce  such  an  effect, 
with  blast  damage  extending  out  to  radii  on  the  order  of  6  km  (and 
fire  to  around  12  km)  [Johnson  and  Larson,  1982],  The  height  of  the 
fuel  zone  is  thus  simply  assumed  constant  and  equal  to  its  outer-belt 
value  in  the  baseline  case.  The  total  height  of  the  combustion  zone 
is  chosen  similarly. 

The  blast-modified  areal  heating  rates  are  not  considered  inde¬ 
pendent  of  radius,  however.  In  each  case,  the  combustibles  of  the 
central  city  and  inner  belt  would  be  spread  radially  to  some  degree 
by  the  blast,  and  piled  up  in  a  debris  field.  Since  some  combustibles 
in  that  zone  may  be  buried  under  layers  of  nonflammable  materials 
(e.g.,  concrete,  brick,  metal),  the  areal  heating  rate  is  not  expected 
to  be  correspondingly  higher,  and  may  in  fact  be  relatively  small. 

We  consider  such  a  situation  in  each  simulation.  It  is  assumed  that 
the  areal  heating  rate  is  zero  at  the  fire  center,  increases  linearly 
with  radius  over  the  debris  zone,  and  equals  its  baseline  outer-belt 
value  for  radii  greater  than  6  km. 
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In  each  model  city  simulation,  the  required  data  set  consists  of 
the  geometric  parameters  R  and  H,  the  heating  rate  parameters  QH  and 
q(r,  v) ,  the  radiation  coefficient  a,  and  the  turbulence  coefficients 
and  K  .  Once  R,  H,  and  QH  are  chosen,  the  parameter  A  in  Eqs.  (1) 
is  determined  along  with  the  radial  velocity  scale  U  from  Eqs.  (2). 

In  all  cases,  the  fire  radius  R  is  taken  to  be  12  km.  However, 
since  the  combustion  zone  height  is  not  constant  in  the  baseline  cases, 
there  are  no  natural  choices  for  H.  The  values  for  H  listed  in  Table  8 

Table  8 

PARAMETERS  IN  MODEL  CITY  SIMULATIONS 


City 


Parameter  W  M  E 


Fundamental  Parameter's 


R 

(km) 

12 

12 

12 

H 

(m)  ^ 

40 

60 

100 

QH 

(kcal/m--sec) 

9.50 

14.25 

23.75 

Auxiliary  Parameters 


U  (m/sec) 

33.6 

33.6 

33.6 

A 

0.347 

0.521 

0.868 

o 

0.264 

0.264 

0.264 

Mr  Kl 

0.125 

0.125 

0.125 

are  arbitrarily  selected  as  characteristic  of  the  baseline  geometry 
in  Figs.  24  through  26.  The  values  of  the  heating  rate  scale  QH  in 
Table  8  are  chosen  likewise. 

Table  8  also  lists  the  corresponding  values  of  U  and  A,  calcu¬ 
lated  using  Eqs.  (2).^  From  Eqs.  (2)  and  (3),  0  depends  on  H,  QH, 
and  the  radiation  mean  free  path  k*  and  depends  on  R,  U,  and 

the  dimensional  eddy  viscosity  2, /p  .  As  in  our  earlier  simulation 

1  a 

j- 

The  value  of  the  dimensionless  perturbation  pressure  scale  6 
is  1.44  x  10~2  f0r  all  cases. 


r 
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(Part  I),  we  take  k*  *  20  m,  2  / p  »  5.04  x  104  m2/sec,  and  M,  =  K, . 

X  a  11 

The  values  we  use  for  a,  M^,  and  [from  Eqs.  (2)  and  (3)]  are  also 
listed  in  Table  8. 

Finally,  we  use  the  heating  rate  distributions  q(r,  y)  defined 
in  Figs.  27  through  29.  The  heating  rate  is  taken  to  be  constant  in 
the  fuel  zone  of  each  baseline  city  region,  and  to  decrease  linearly 
with  altitude  over  the  remainder  of  the  combustion  zone,  approaching 
zero  at  the  top  of  that  zone.  Subject  to  the  choices  of  QH  in  Table  8, 
the  fuel  zone  values  of  q(r,  y)  are  selected  so  that  the  resulting 
areal  heating  rates 


QH  q(r ,  y)  dv 

are  equivalent  to  those  in  Figs.  24  through  26.  The  blast-modified 
forms  of  q(r,  y)  are  developed  likewise. 

SIMULATION  RESULTS 

The  results  of  the  model  city  simulations  are  summarized  by 
Table  9,  which  compares  predicted  baseline  and  blast-modified  velocity, 
temperature,  and  pressure  maxima.  As  expected,  the  predictions  for 
the  baseline  simulations  are  all  larger  than  those  for  the  blast- 
modified  simulations.  The  differences  are  great  for  the  temperature, 
pressure,  and  vertical  velocity  maxima,  but  small  for  the  radial 
velocity  maxima  (10  percent  or  less).  Therefore,  the  winds  and  wind 
damage  resulting  from  nuclear-weapon-ignited  fires  may  be  relatively 
insensitive  to  the  degree  of  blast  disruption  of  the  fuel  bed. 

In  any  case,  the  predictions  in  Table  9  indicate  that  the  winds 
generated  by  a  large  urban  fire  will  in  themselves  constitute  a  major 
threat.  Although  most  of  the  velocities  in  the  table  are  less  than 
hurricane  force  (more  than  30  m/sec),  it  should  be  noted  that  those 
values  represent  means.  Near  street  level,  where  fire  winds  will  be 
channeled  between  buildings,  wind  speeds  of  hurricane  force  should 
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iii.  p  i  wm. 


Table  9 

VELOCITY,  TEMPERATURE,  AND  PERTURBATION  PRESSURE 
MAXIMA  IN  MODEL  CITY  SIMULATIONS 


Simulation 

City 

W  M 

E 

Maximum  Radial  Velocity,  u^tox 

Baseline 

20.2  26.3 

39.0 

Blast-Modified 

17.9  23.9 

28.5 

Maximum  Temperature ,  Tmax  (°K) 

Baseline 

577  619 

Blast-Modified 

455  485 

510 

Maximum  Perturbation  Pressure,  Pma~  (psi) 

Baseline 

0.056  0.113 

0.271 

Blast-Modified 

0.011  0.044 

0.076 

Maximum  Vertical  Velocity,  vmax 

(m/sec) 

Baseline 

0.89  3.12 

12.48 

Blast-Modified 

0.37  1.56 

4.32 

be  considered  typical.  The  winds  may  be  even  greater  than  those 
encountered  in  the  Hamburg  firestorm  of  1943  (see  Part  I). 

Also  as  expected,  the  velocity,  temperature,  and  pressure  pre¬ 
dictions  in  Table  9  are  all  greatest  for  city  E  (the  tallest  and 
densest)  and  least  for  city  W  (the  shortest  and  sparsest).  For  a 
given  fire,  therefore,  the  threat  will  be  most  severe  for  the  most 
congested  cities.  In  general,  however,  the  shorter  cities  sprawl 
out  over  greater  areas  than  do  taller  ones  of  comparable  population, 
and  are  thus  capable  of  supporting  more  widespread  fires.  Multiple 
weapon  bursts  can  greatly  increase  the  fire  severity  in  such  cities. 

In  the  remainder  of  this  subsection,  we  describe  the  results 
of  the  model  city  simulations  in  more  detail.  The  vector  plot 
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in  Fig.  30  is  representative  of  the  flow  fields  obtained  in  the  sim¬ 
ulations.  The  flow  in  each  case  is  everywhere  inward  (toward  the 
symmetry  axis)  and  upward,  with  the  horizontal  velocity  typically 
maximum  on  the  symmetry  axis  in  the  upper  portion  of  the  turning 
region.  Figures  31  through  34  present,  respectively,  the  temperature 
and  pressure  contours  and  radial  and  vertical  velocity  profiles  for 
baseline  and  blast-modified  simulations  for  city  W.  Similar  sets  of 
plots  for  cities  M  and  E  are  presented  in  Figs.  35  through  38,  and 
39  through  42,  respectively. 

As  expected  (cf.  Fig.  12),  the  baseline  and  blast-modified  tem¬ 
perature  contours  in  Fig.  31  are  quite  different.  In  the  baseline 
case,  the  heat  release  is  greatest  at  the  center  of  the  fire  (see 
Fig.  27)  and  the  temperature  is  maximum  on  the  symmetry  axis.  At 
any  fixed  altitude,  the  variation  of  temperature  with  radius  is 
roughly  Gaussian.  However,  in  the  blast-modified  case,  the  tempera¬ 
ture  maximum  does  not  occur  on  the  symmetry  axis.  Rather,  since  the 
heat  release  is  greatly  diminished  and  actually  least  in  the  fire 
center  (Fig.  27),  the  temperature  is  maximum  about  halfway  between 
the  center  and  the  periphery.  Additionally,  the  temperature  is  uni¬ 
formly  reduced  in  the  blast-modified  case. 

Figures  32  through  34  show  little  qualitative  difference  between 
the  baseline  and  blast-modified  pressures  and  velocities.  The  con¬ 
tour  plots  in  Fig.  32  exhibit  the  same  behavior,  the  pressure  drop 
decreasing  with  an  increase  in  either  radius  or  altitude.  Likewise, 
the  profile  plots  in  both  Figs.  33  and  34  are  qualitatively  the  same. 
The  radial  velocity  is  maximum  at  ground  level  on  the  periphery  and 
the  vertical  velocity  is  maximum  on  the  symmetry  axis  in  the  upper 
part  of  the  turning  region.  At  any  specified  altitude,  the  variation 
of  vertical  velocity  with  radius  is  again  roughly  Gaussian.  Of  course, 
as  expected,  certain  quantitative  differences  are  evident:  pressure 
drops  and  velocities  are  all  less  in  the  blast-modified  case  than  in 
the  baseline  case.  However,  the  radial  velocity  differences  are  rela¬ 
tively  small. 

The  comparisons  in  both  Figs.  35  through  38  and  Figs.  39  through 
42  are  similar  to  those  in  Figs.  31  through  34.  The  primary  difference 
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Fig.  32--Baseline  and  blast-modified  pressure  contours 

for  city  W 
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Fig.  33--Baseline  and  blast-modified  radial  velocity  profiles 

for  city  W 


35--Base1ine  and  blast-modified  temperature  contours 

for  city  M 
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Pig.  36--Baseline  and  blast-modified  pressure  contours 

for  city  M 
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Fig.  40--Base1 ine  and  blast-modified  pressure  contours 

for  city  E 
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Fig.  41--8aseline  and  blast-modified  radial  velocity  profiles 

for  city  E 


in  the  results  presented  in  the  corresponding  plots  is  that  temper¬ 
atures,  pressure  drops,  and  velocities  all  increase  with  total  load¬ 
ing,  and  are  least  for  city  W  and  greatest  for  city  E.  Further,  as 
shown  in  Fig.  41,  the  baseline  radial  velocity  for  city  E  is  not 
maximum  at  the  fire  periphery,  as  in  the  other  cases  (cf.  Figs.  33 
and  37),  but  rather  about  halfway  between  the  periphery  and  the  center. 
That  difference  presumably  reflects  the  greater  radial  extent  of  the 
taller  central-city  and  inner-belt  regions  in  this  case  than  in  the 
others  (cf.  Figs.  27  through  29). 


IV.  SAMPLE  TIME-DEPENDENT  SIMULATION 


This  section  extends  the  quasi-steady  turning-region  analysis 
of  Part  I  to  obtain  time-dependent  predictions.  The  multiple- fuel- 
bed  Flambeau  fire  [Countryman,  1969;  Palmer,  1981]  treated  in  Part  I 
is  used  as  a  sample  case. 

Figure  43  presents  a  typical  time-history  of  the  areal  heat  re¬ 
lease  in  the  largest  Flambeau  fires.  During  the  period  encompassing 
most  of  the  burning  (5  to  55  min  after  the  fire  starts),  variations 
in  the  mean  heat  release  rate  occur  on  a  time  scale  of  many  minutes. 
At  any  given  point  during  that  period,  a  steady-state  analysis  like 
that  in  Part  I  (Sec.  IV)  is  justified  as  long  as  dynamic  transients 
occur  on  a  much  shorter  time  scale  of,  say,  a  minute  or  less.  As 
we  now  show,  that  criterion  is  satisfied  at  all  times  between  5  and 
45  min.  A  time-dependent  simulation  is  constructed  by  performing 
steady-state  calculations  at  several  points  in  that  interval. 


Time  (min  after  fire  starts) 

Sawea:  Adaptad  from  Fig.  3  of  Faknar  (1981).  Raprintad  with  parmMen  from  A tmotphtric 
Environmtnt,  Vel.  IS,  Thamat  Y.  Palmar,  "Largo  Fira  Wind*.  Gaaaa  and  Smafca,"  1981,  Fargamen 
Praaa.  Ltd. 


Fig.  43--Areal  heat  release  time-history  for  Flambeau  fire  760-12 


In  the  steady-state  Flambeau  simulation,  we  used  R  =  250  m, 

2 

H  =  20  m,  QH  ~  50  kcal/m  -sec,  and  hence  had  U  ~  8  m/sec  [see  Eqs. 
(2)].  For  that  case,  the  time  scale  R/U  on  which  dynamic  transients 
occurred  was  approximately  30  sec,  satisfying  the  criterion  for  the 
flow  to  be  quasi-steady.  In  the  new,  time-dependent  analysis,  we 
similarly  use 


R  =  250  m  , 


H  =  20  m  , 


(10) 


but  vary  QH  with  time  as  shown  in  Fig.  A3.  At  any  given  time,  the 

9 

criterion  for  the  flow  to  be  quasi-steady  is  that  QH  >,  25  kcal/m“-sec . 
From  Fig.  43,  it  is  certainly  satisfied  for  periods  between  5  and 
45  min,  and  possibly  up  to  55  min  as  well. 

In  the  following  simulation,  we  calculate  quasi-steady  flows  at 
points  5  min  apart.  The  specific  values  of  QH  used  at  those  points 
are  summarized  in  Table  10.  As  in  the  steady-state  simulation  of 
Part  I,  we  also  use 


/1.6 

for 

0  S  y  S  0.25 

q(r,  y)  =  < 1.6 

(i  (1  -  y>) 

for 

0.25  Sy  £  1.0  , 

(o 

for 

y  a  1.0 

*-l 

k  =20  m  , 

Z,  k  2 

—  =  -  =  8410  m  /sec  (11) 


at  each  point.  Thus,  from  Eqs,  (2)  and  (3), 


U  =  8.41  x  F  m/sec  , 


2.77 
_2  ’ 


0.022 

F 


4.0 

al  F  * 


where 


57  kcal/m  -sec 

and,  from  p.  6,  6  =  9.04  *  10  ^  x  F^. 

Figures  44  through  47  summarize  the  results  of  the  time-dependent 
Flambeau  simulation.  The  predicted  variations  with  time  of  the  tem¬ 
perature,  pressure,  and  velocity  fields  are  shown  for  the  heat  release 
described  in  Fig.  43.  The  solid  curves  represent  the  results  of  the 
quasi-steady  analysis,  which  is  considered  appropriate  for  the  period 
approximately  5  to  55  min  after  the  fire  starts. 

Table  10 

HEAT  RELEASE  SCALES  USED  IN  TIME- 
DEPENDENT  FLAMBEAU  SIMULATION 


Time  (min) 

QH  (kcal/m^-sec) 

5 

134 

10 

84 

15 

53 

20 

38 

25 

31 

30 

24 

35 

20 

40 

18 

45 

18 

50 

14 

55 

10 
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Fig.  46--Time-hi story  of  induced  peripheral  fire  wind  (radial  velocity) 

for  sample  Flambeau  fire 


Time  (min) 

Fig.  47--Time-bistory  of  emerging  column  flow  (vertical  velocity  at  top 
of  turning  region)  for  sample  Flambeau  fire 
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During  the  first  5  min  of  burning,  when  the  individual  fires 
are  coalescing  to  form  a  large-fire  system,  the  heat  release  rate 
(Fig.  43)  is  changing  on  a  time  scale  that  is  fairly  fast  compared 
with  that  for  convection.  The  same  is  also  true  during  the  final 
phase  of  fire  activity,  i.e.,  the  period  later  than  1  hr  or  so  when 
the  large  fire  begins  to  disintegrate  and  burn  out.  For  those  two 
periods,  the  quasi-steady  analysis  is  thus  inappropriate,  and  the 
dashed  curves  in  Figs.  44  through  47  illustrate  the  expected  results 
of  a  more  rigorous  analysis. 

Figure  44  describes  the  predicted  variation  of  the  temperature 
field  with  time  in  the  burning  zone.  The  contours  shown  in  Fig.  48 
are  representative  of  the  temperature  field  at  any  given  time.  Maximum 
temperatures  occur  in  the  center  of  the  fire  just  above  the  top  of 
the  fuel  zone  (height  «  4  m) .  Minimum  temperatures  occur  on  the  per¬ 
iphery  at  the  top  of  the  burning  region.  As  expected,  the  time- 
histories  of  those  temperatures  resemble  the  heat  release  time-history 
shown  in  Fig.  43.  Temperatures  are  everywhere  greatest  about  5  min 
after  the  fire  starts,  when  the  heat  release  is  greatest,  and  then 
decrease  with  time. 

The  corresponding  pressure  and  velocity  time-histories  shown  in 
Figs.  45  through  47  follow  the  same  pattern.  As  illustrated  in  Figs. 

49  and  50,  pressure  gradients  and  radial  velocities  at  any  given  time 
are  maximum  at  the  fire  periphery.  At  all  times,  they  are  greatest 
at  ground  level  and  decrease  with  altitude.  Finally,  Fig.  47  char¬ 
acterizes  the  variation  with  time  of  the  vertical  velocity  emerging 
at  the  top  of  the  computational  region  (10  flame  heights,  200  m) .  As 
illustrated  in  Fig.  51,  velocity  profiles  at  that  height  are  almost 
top-hat  at  all  times.  Maximum  and  minimum  velocities,  which  occur 
above  the  fire  center  and  on  the  periphery,  respectively,  are  thus 
always  quite  similar. 
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V.  DISCUSSION 


In  presenting  calculations  of  the  urban  fire  environment  that 
may  result  from  a  nuclear  weapon  burst,  our  aim  has  been  to  describe 
the  macroscopic  features  of  a  large  area  fire  through  analysis  of  the 
interactions  that  produce  high  temperatures  and  fire  winds.  Our 
approach  has  been  to  develop  a  model  based  on  first  principles.  In 
contrast  to  asymptotic  theories  that  extrapolate  the  properties  of 
a  free-convection  plume,  we  consider  the  fire  region  directly.  De¬ 
spite  our  use  of  several  simplifications  and  empiricisms,  we  expect 
that  the  theory  yields  a  fine  resolution  of  the  hydrothermodynamics 
of  a  large  area  fire. 

The  analysis  considers  a  quasi-steady  axisymmetric  fire  in  an 
urban  area.  A  spatially  dependent  volume-heat-release  function  is 
used  to  describe  the  energy  input  from  fires  in  a  blast-disturbed 
urban  fuel  bed.  Though  the  function  can  be  defined  to  reflect  vary¬ 
ing  fuel  loads  and  distributions  of  burning,  the  results  presented 
are  based  on  sectionally  uniform  heating  rates. 

The  formulation  of  the  boundary  value  problem  is  unique  in  that 
analytic  jump  conditions  are  derived  for  the  fire-column  boundary. 

Such  a  formulation  obviates  the  need  for  extensive  far-field  calcula¬ 
tions  and  allows  demonstration  of  the  dependence  of  induced  fire  winds 
on  the  fire  area  and  heat  release  rate.  Improvements  to  the  model 
could  include  the  use  of  multiparameter  models  to  describe  the  struc¬ 
ture  of  local  turbulence  and  hot-gas/smoke  radiation  (in  place  of  the 
simple  one-parameter  formulations  used  here)  as  well  as  a  heat  re¬ 
lease  function  dependent  on  the  local  thermodynamic  and  flow  conditions. 

The  effects  of  the  model  parameters  were  explored  in  a  sensitivity 
study.  Varying  fire  areas,  heat  distributions  and  release  rates,  mix¬ 
ing  coefficients,  and  radiation  lengths  were  investigated.  Qualita¬ 
tively  similar  flows  were  observed  for  all  parameter  ranges,  despite 
significant  quantitative  differences.  The  fire-induced  velocities 
increased  with  both  heat  release  and  fire  area.  For  small  values,  the 
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velocity  increased  almost  linearly.  For  larger  values  of  either  fire 
area  or  heat  release  rate,  the  induced  velocities  increased  more 
slowly.  Presumably,  the  magnitudes  are  limited  by  compressibility 
effects  through  a  slowing  production  of  buoyancy. 

Variations  in  the  degree  of  turbulent  mixing  produced  substan¬ 
tial  changes  in  the  velocity  fields.  As  expected,  large  mixing  coef¬ 
ficients  tend  to  reduce  the  velocities  at  the  fire  periphery;  smaller 
values  produce  higher  velocities.  In  general,  the  mixing  coefficients 
should  reflect  both  the  apparent  surface  roughness  (urban  structure) 
and  the  fire-generated  turbulence.  The  present  analysis  used  constant 
eddy  dif fusivities . 

Fuel  distributions  and  heat  release  rates  were  developed  to  ex¬ 
plore  the  effects  of  varying  city  constructions  and  fuel  densities  on 
the  fire  environment  in  three  model  urban  areas.  Though  the  models 
provide  only  a  simple  representation  of  a  city  cross  section,  the 
results  indicate  the  importance  of  modeling  a  particular  fuel  bed. 

More  refined  distributions  based  on  city  surveys  could  be  used. 

A  fire  radius  of  12  km  is  assumed  for  each  model  urban  area, 
since  megaton-yield  explosions  would  cause  area  fires  of  that  magni¬ 
tude.  Each  model  area  is  characterized  by  a  high-density  city  center, 
a  lower  density  annulus  representing  a  mixed  residential/ industrial 
construction,  and  an  outer  belt  of  low-density  residential  tracts. 

The  heat  release  function  reflected  the  sectionally  variable  fuel 
densities  and  loadings. 

Computed  result.'  for  each  fire  illustrate  the  dependence  of  the 
fire  winds  and  temperature  levels  on  the  urban  geography.  The  most 
severe  fires  occur  in  the  higher  density  cities,  though  even  the  lower 
density  constructions  support  high  temperatures  and  velocities. 
Application  of  those  results  to  definition  of  shelter  hardness  (ther¬ 
mal)  would  imply  different  criteria  for  each  city. 

An  additional  set  of  numerical  experiments  considered  each  model 
area  with  a  reduced  level  of  burning  in  the  center  region.  A  slow- 
burning  debris  field  created  by  extensive  blast  damage  may  foster  such 
a  heat  release  distribution.  In  general,  peak  velocities  were  only 
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slightly  reduced,  though  temperature  changes  were  substantial.  Tem¬ 
perature  levels  in  the  central  region  were  markedly  reduced  from  those 
in  the  full-burning  outer  annulus.  It  is  reasonable  to  expect  that 
the  fire  winds  would  fan  the  slow-burning  central  area,  creating  a 
more  severe  environment . 

The  analysis  and  results  presented  here  describe  the  environment 
of  a  large  area  fire  at  fixed  times.  Past  experience  indicates  that 
such  fires  develop,  peak,  and  ebb  over  a  period  of  several  hours. 

The  methods  developed  here  can  be  used  successively  to  estimate  the 
environment  at  different  times  in  the  fire  evolution.  The  high  ve¬ 
locities  and  temperatures  predicted  would  encourage  a  rapid  internal 
fire  spread,  which  may  considerably  hasten  burnout  of  the  city.  Im¬ 
proved  estimates  would  result  from  use  of  a  responsive  heat  release 
function  that  depends  on  the  local  gas  dynamic  state.  A  time-dependent 
calculation  that  includes  the  high-resolution  analysis  presented  here 
could  provide  a  complete  map  of  the  evolution  and  ebb  of  urban  fires. 
Synthesis  of  the  present  analysis  with  a  large  hydro-code  may  provide 
a  numerically  efficient  model  for  calculating  the  m esoscale  flow  field 
of  a  large  area  fire. 
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Appendix 

PREDICTION-ALGORITHM  DOCUMENTATION 


This  appendix  documents  the  computer  algorithm  used  to  generate 
the  predictions  and  other  results  presented  in  Secs.  II  through  IV. 
The  algorithm  yields  numerical  solutions  to  the  turning-region 
boundary  value  problem  derived  in  Part  I. 

The  turning-region  problem  is  defined  by  the  following  balances: 


^  (rpu)  +  (rpv)  =  0  , 
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f  +  AP  =  0  . 
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The  flow  chart  In  Fig.  A.l  summarizes  the  solution  algorithm.  As 
indicated  in  the  figure  (and  discussed  in  Part  I,  pp.  16-17),  the 
algorithm  uses  an  iterative  "shooting"  method  [Keller,  1968]  to  solve 
the  overall  boundary  value  problem.  Each  shoot  employs  a  nonlinear 
Crank-Nicolson  finite  difference  scheme. 

INPUTS  AM)  OUTPUTS 

The  following  is  an  ordered  list  of  the  FORTRAN  inputs  required 
for  automated  algorithm  usage: 


a.  System  parameters: 

A  A 

SIGMA  a 
XNUIN  M1 
XKIN  K 

b.  Mesh  descriptors: 


YMAX 

M 

N 


y 

max 

Number  of  finite  difference  cells  in  the  radial 
direction. 

Number  of  finite  difference  cells  in  the  vertical 
direction. 1 


c.  Iteration  parameters: 


ETLINE  Maximum  absolute  error  allowed  in  solving  dis¬ 
cretized  Crank-Nicolson  equations  on  any  line  of 
constant  y. 

ETOLP  Maximum  absolute  error  allowed  in  P  on  y  =  v 

■  max 

MAX1TL  Maximum  number  of  Newton  iterations  allowed  in 

solving  discretized  equations  on  any  line  of  con¬ 
stant  y. 

MAXITP  Maximum  number  of  complete-shoot  iterations  allowed 
in  solving  the  complete  boundary  value  problem. 


the  present  coding,  cell  widths  are  taken  to  be  uniform,  the 

common  width  being  denoted  by  Ar.  For  0  s  y  s  5,  cell  heights  Ay  are 

also  taken  to  be  Ar;  for  5  s  y  s  y  ,  thev  are  taken  to  be  2  Ar. 

max 


i 


V  sufficiently  near 


Generara  new  approximation  of  ?(r»  0)  o*ng  simple  iteration 

P(r.  0)  _  -  Pf*.  0)  -  u  ( ?(r,  y  )  ♦  Ay  )  , 

nev  old  \  max  max  / 


being  rtiixftion  ftetor 


Fig.  A.l--Flow  chart  of  prediction  algorithm  for  turning-rec 

boundary  value  problem 


d.  Program  control  parameters: 


(Qq)__  .  maximum  value  of  Qq  in  present  coding. 
Parameters  available  for  either  system  descrip¬ 
tion  or  program  control.  (These  parameters 
are  currently  unused,  so  they  are  set  equal 
to  zero.) 

Parameter  that  controls  the  initial  choices 
of  first-shoot  values  along  y  =  0. 

IGUESS  *  1:  Discretized  versions  of  P(r,  0), 
T(r,  0),  and  u(r,  0)  must  be  sup¬ 
plied  as  inputs. 

IGUESS  =  2:  "Standard"  initial  choices  are 
made  automatically. 

IGUESS  =  3:  Discretized  versions  of  P(r,  0), 
T(r,  0),  and  u(r,  0)  from  the 
last  algorithm  run  are  used  auto¬ 
matically. 

The  "standard"  initial  guesses  in  the  present  coding  represent  the 
ground-level  pressure,  temperature,  and  radial  velocity  profiles 
predicted  in  Part  I  (pp.  13-16)  for  the  case  of  a  weakly  heated  flow: 

°>  -  -  «*w(£)  • 

T(r,  0)  -  1.0  +  (*,>„,„(£)  , 

u(t, 0)  -  -  (^)Mx  (»sq;) r  •  w-3’ 

Finally,  the  functional  form  of  q(r,  y)  must  be  prescribed. 
Currently,  it  is  given  in  the  form  used  in  our  analysis: 
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((Qq)max 

q(r,  y)  =  <(Qq)max(f  d  -  X>) 


for  0  s  y  s  0.25 

i 

for  0.25  £  y  £  1.0 

1 


(A. 4) 


Changes  require  reprogramming  in  subroutine  BCFUNC.  As  a  user  option, 
the  dissipation  coefficients  and  may  also  be  made  to  depend  on 
r  and  y  through  reprogramming,  also  in  BCFUNC.  Currently,  those  co¬ 
efficients  are  simply  taken  as  constants. 

A  sample  input  data  stream  is 

A  =  0.217  ,  SIGMA  =  0.110  ,  XNUIN  =  XKIN  =  0.200  ; 

YMAX  =  10  ,  M  »  40  ,  N  =  300  ; 

ETLINE  =  0.0003  ,  ETOLP  -  0.001  ,  MAXITL  «  30  ,  MAXITP  =  9  ; 

ALPHA  =  1.6  ,  BETA  =  0.0  ,  BLOSS  =  0.0  ,  IGUESS  -  3  '.  (A. 5) 

The  corresponding  output  is  listed  in  Table  A.l.  All  output  data  are 
automatically  dumped  onto  a  disk  file  and  used  in  constructing  the 
vector,  contour,  and  profile  plots  presented  in  Secs.  II  through  IV. 

PROGRAM  /LOW 

Figure  A.l  outlines  the  overall  flow  of  the  computer  algorithm. 
Inputs,  outputs,  and  the  basic  shooting  iteration  are  controlled  by 
the  main  program.  Figure  A. 2  presents  a  flow  chart  of  that  program. 

In  subroutine  BCFUNC,  functional  forms  that  may  be  varied  [e.g., 
q(r,  y)]  are  specified;  and  in  subroutine  SPRINT,  final  output  data 
are  printed.  Subroutine  MSWEEP  performs  most  of  the  computations — 
l.e.,  those  providing  a  one-shoot  finite  difference  solution  to 
Eqs.  (A.l)  for  given  data  along  y  =  0,  r  =  0,  and  r  =  1.  The  flow 
of  that  subroutine  is  shown  in  Figs.  A. 4  and  A. 5  below  (pp.  80-81); 
the  underlying  numerical  analysis  is  detailed  in  the  accompanying 


SAMPLE  CODE  OUTPUT 


Table  A.l — Continued 
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text.  Complete  FORTRAN  listings  of  the  main  program  and  all  sub¬ 
routines  are  presented  at  the  end  of  this  appendix. 

FINITE  DIFFERENCE  ALGORITHM  AND  SUBROUTINE  MSWEEP 

The  finite  difference  scheme  underlying  the  computer  algorithm 
is  based  on  the  rectangular  grid  and  "stencils”  in  Fig.  A. 3.  The 
grid  has  uniform  cell  widths  (Ar  *  1/M)  and  piecewise  uniform  heights 

(Ay  =  Ar  for  0  £  y  £  5,  y=2Ar  for  y  >  5).  Mesh  points  are  thus 

located  at  (r^,  y^)  points  with 

r±  *  (i  -  1)  Ar  ,  1  £  i  £  M  +  1  , 

(  (j  -  1)  Ar  ,  1  £  j  £  5M  +  1  ) 

=  j  }  •  (A. 6) 

J  f  5  +  (j  -  SM  -  1) (2  Ar)  ,  5M  +  2  £  j  £  N  +  1  j 

Dependent  variables  evaluated  at  those  points  are  given  corresponding 

(i,  j)  suffices.  For  example,  r  *  0  (and  not  Ar) ,  y  =  y 

i-  N '  i.  ma  x 

and  =  u(0,  y^^) .  The  accurate  approximation  of 

{ui,j’  vi,j’  Pi,j’  Ti,j*  Pi,j}l£i£M+l,l£j£N+l  constitutes  a  numerical 
solution  of  the  problem  posed  by  Eqs.  (A.l)  and  (A. 2).  The  stencils 
in  Fig.  A. 3,  explained  below  (pp.  82-83),  indicate  the  types  of  dif¬ 
ferencing  used  in  the  algorithm. 

Along  the  y  *  0  line  (where  j  -  1),  v  =  0.  As  shown  in  Fig.  A. 2, 

u  i  1 

a  guess  at  P  along  that  line  [actually,  the  (P.,}‘  7  values]  is  also 

il  1 

made  before  the  MSWEEP  subroutine  is  called.  As  indicated  in  the 
first  large  block  of  Fig.  A.l,  Eqs.  (A.l)  and  (A. 2)  then  reduce 
along  the  y  *  0  line  to  the  following  ordinary-differential-equation 
boundary  value  problem  for  u  and  T  alone  (with  p  =  1/T) : 


* 
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and  the  components  of  F  are  the  finite  difference  equations  obtained 
from  the  discretization  of  Eqs.  (A. 7).  The  Crank-Nicolson  solution 
of  Eqs.  (A.l) — subject  to  Eqs.  (A. 2b)  and  (A. 2d),  and  the  solution 
of  Eqs.  (A. 7)  [with  P(r,  0)  given] — also  involves  the  solution  of 
Eq.  (A. 8)  for  successive  j,  j  *  2,  3,  4,  ...»  N,  N  +  1.  In  that  case, 
the  components  of  F  are  slightly  different  from  those  for  j  *  1.  For 
each  fixed  j (a2) ,  x  is  redefined  [cf.  Eq.  (A. 9)]  as 


The  overall  flow  in  subrout^.c  MSWEEP  is  thus  as  outlined  in  Fig.  A. 4, 
Newton's  method  (see  Fig.  A. 5)  being  used  in  all  cases  to  solve  the 
component  equations  of  Eq.  (A. 8)  because  they  are  nonlinear. 

To  fully  describe  subroutine  MSWEEP,  we  must  define  the  follow¬ 
ing  processes: 


1.  The  specific  finite-difference  schemes  used  to  reduce  Eqs. 

(A. 7)  to  Eq.  (A. 8)  for  j  =  1  (y  *  0),  and  Eqs.  (A.l),  (A. 2b), 
and  (A. 2d)  to  Eq.  (A. 8)  for  j  i  2  (y  >  0) . 

2.  The  corresponding  specifications  of  J  (as  9F/3x;  see  Fig. 

A.  5) . 

3.  The  specific  Gaussian  elimination  scheme  used  to  compute 
5  (as  -J-1F ;  see  Fig.  A. 5). 


r 


Fig.  A. 4— How  chart  of  subroutine  MSWEEP 
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Fig.  A. 5— Flow  chart  of  Newton  steps  used  in  subroutine  MSWEEP 
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For  j  -  1,  che  usual  discretization  for  second-order  ordinary  dif¬ 
ferential  equations  is  used  to  reduce  Eqs.  (A. 7)  to  Eq.  (A. 8).  That 
is,  for  any  dependent  variable  w,  interior  derivatives  are  replaced 
by  the  centered  difference  approximations 


w  —  2  v  +  w 

i+l.l  i,l  i-l,l 


(Ar)' 


/  3w\  ^  Wi+l.l  ~  wi-l,l 

\3r4,i  2  Ar 


(A. 11) 


In  addition,  boundary  derivatives  are  taken  as,  for  example. 


(A. 12) 


Thus,  as  indicated  by  the  corresponding  discretization  stencil  at 
the  bottom  of  Fig.  A. 3,  no  more  than  three  values  of  u  and  T  are 
involved  in  each  component  equation  in  Eq.  (A. 8)  [see  Eq.  (A. 9)] — 
and  those  values  are  at  successive  r^.  For  each  j  a  2,  the  Crank- 
Nicoison  scheme  to  be  defined  by  Eqs.  (A. 13)  and  (A. 14)  is  used  to 
reduce  Eqs.  (A.l),  (A. 2b),  and  (A. 2d)  to  Eq.  (A. 8).  Centered  dif¬ 
ferences  are  also  used,  with  interior  and  boundary  approximation 
stencils  as  shown  in  Fig.  A. 3,  interior  approximations  being 


i 
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and  the  boundary  formulas  being  of  the  form 


(A. 13) 


As  in  linear  problems  [Isaacson  and  Keller,  1966],  the  resulting 
nonlinear  scheme  has  proven  to  be  convergent,  stable,  and  accurate 
in  providing  solutions  to  the  turning-region  boundary  value  problem. 

In  reducing  either  Eqs.  (A. 7)  or  Eqs.  (A.l),  (A. 2b),  and  (A. 2d) 
to  Eq.  (A. 8)  by  means  of  Eqs.  (A. 11)  through  (A. 14),  it  proves  con¬ 
venient — both  notationally  and  computationally — to  begin  by  decom¬ 
posing  x  (for  all  j)  into 


f  * 

/-  ' 

Vi 

T,  . 

2,J 

U3,i 

T-  • 

3»  J 

• 

— ► 

• 

u  * 

• 

and  T  = 

• 

Vi 

T 

M,j 

“M+l,j 

T 

M+l,j 

and  arranging  the  corresponding  F(x)  as 


(A. 15) 
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vhere  the  and  represent  the  components  that  arise  primarily 
from  the  discretization  of  the  horizontal  momentum  and  energy  equa¬ 
tions,  respectively,  about  the  grid  point  (r^  y^ ) .  The  specific 
forms  of  the  and  are  derived  shortly.  As  mentioned  above 
(p.  82),  centered  differencing  results  in  each  G^  and  involving 

at  most  six  of  the  2M  u,  .  and  T.  ..  The  Jacobian  matrix  derivatives 

i.J  i,J 


with  respect  to  u  and  ^  thus  have  the  following  sparse,  banded 
forms : 


kwgrJlXt*m* 


From  Eqs.  (A. 10)  and  (A. 15)  through  (A. 18),  the  J  -  3F/3x  Jacobian 
is  then  easily  computed  from  the  nonzero  elements  in  the  Jacobians 
of  Eqs.  (A. 13)  as 


=  u  *32  it: 


'31..  ’32 


">;.(>  V;.i  V-2.2  *1-2,2 


vn-2,;  5m-:.i  '*-2.2  :,,-2.2 


N-i.i)  "v-i.o  N-i.i  *1-1.1  *1- 


:«-;.o  Vi.i  'M-i.i 
’  V,  S.n 


DIFFERENCE  EQUATIONS  AND  JACOBIAN  ELEMENTS 

Specific  forms  for  G  and  H,  and  hence  for  F  and  the  various  a, 
S,  y,  and  5  coefficients  in  Eqs.  (A. 18)  and  (A. 19),  are  derived  as 
follows.  Subject  to  Eqs.  (A. 11),  the  first  two  lines  of  Eqs.  (A. 7) 
are  approximated  at  r  ■  r^,  2  S  i  £  M,  by 
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H  »  h,  (T  -  T  )  -  Um+1>1 

^i+1  Ar  'TH-1,1  M,1  T_.  _  ViM+l,l 

cirri ,  1 


(A. 


From  Eqs.  (A. 23)  through  (A. 25),  the  Jacobian  elements  in  Eqs. 
(A. 18)  and  (A. 19)  are  then  found  to  be 


25) 
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for  3  Si  sX;  and 


M, 


“M,0  Ar 


“».i '  ^  '  2(vli.i)(1  *  Vi.i)  ’ 


^M,0  =  °  ’ 


BM, 


1  \T, 


M+i,iy 


M+l 


a) 


YM,0  “  1  +  T, 


M+1,1 


YM, 1  "  °  * 


M,0  '  Ar  ’ 


M,1  Ar 


UM+1,1 

T2 

M+1,1 


(A. 


For  j  &  2,  Che  derivation  of  relevant  forms  for  G  and  H  begins 
with  a  preliminary  rewriting  of  Eqs.  (A.l)  in  the  following  equiv¬ 
alent  form: 


r  £  (rPu)  +  Jv  (pv)  =  0  * 


7  -h  (r°u2)  +  fr  +  ^v  (PUV) 


M  j  LR  +  I  _  JL 
Hsr2  r2 


28) 


f  +  ap  -  o  , 
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+  q(r,  y)  -  a(T4  -  l)  , 


pT  «  1  . 


(A. 29) 


Subject  to  Eqs.  (A. 13),  this  set  is  approximated  at  each  interior 
grid  point  (r^  y  ),  2  s  i  s  H,  of  line  j,  2  S  j  s  N  +  1  by  the 
finite  difference  system 
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Ay  (Pi.j  -  '  '*((^1,1  +(l)lfJ-i)  ' 


(A. 30c) 


1  /  i+1 


%  i-1  / 

, ) - T —  a  +  u. 


4  Ar  ^  v“i+l,j  '  “i+l.j-l'  r±  r  i-l,j-l' 


+  Ay  Cvi,J  - 


T .  ,  .  -  2T.  .  +  T  ,  .  +  T.^.  .  -  2T.  +  T.  n  .  . 

|  *±Li  1.1  1-1. -i  i+i. -i-i  i.j-i  i-i.J-i 

2(Ar)2 


(T  -  T  +  T  -  T 

Li+1, i  "i-l.i  i+l.i-1  i-1,1-11 

4  Ar  , 

*  I  (’<V  S':1  +  i(Ti-  +  °^(  1,J  21,X'J)  -  ■ 


(A. 30d) 


p  having  been  eliminated  for  convenience.  This  system  is  simplified 

by  solving  Eqs.  (A. 30a)  and  (A. 30c)  for  v.  .  and  P.  .in  terms  of  the 

1  *  J  1  f  J 

various  u  and  T  on  line  j  and  variables  on  line  j  -  1  as 
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where 


li  ■  2  L  '"l+l.J-l  -  2ui,3-l  +  “t-l.j-l’  -  4  (  T1  (t)1+1  ^ 

ri-l  (a2\  _“l,  %  /M 

ri  '  T  4-1  j-1  ri  Ui+l.j-l  ~  ui-l, j-1  +  2yr2 J  i,j-l 


Ar 


+  2(P  -  p  )  .  »  / — 

1  i+i.j-1  i-i,j-r  A  2  \  t . 


i+1, j-1  i-1, j-1 


Ar  / uv\ 

Ay  \ T  4, j-i 


K, 


‘i  2  Ar  '*i+l,j-l  “  2Ti,j~l  +  Ti-l,J-l)  “  4  ((  I*/  Ui+l,j-l 


H..  -  t  TT  (T . 


(T 


-  T  )  I  +  Al  v 

L  -I  1  1  "l  '  M  '  A  __  V4  ■ 


rt  /  i-1, j-1  rt  i+1, j-1  li-l,j-r/  Ay  vi,j-l 


2  (q(ri’  yj>  +  q(V  yj-l>)  Ar 


(A. 33) 


As  outlined  in  Fig.  A. 4,  a  recursive  Crank-Nicolson  scheme  is 
employed  in  subroutine  MSWEEP  to  generate  a  numerical  solution  to 
Eqs.  (A.l),  (A. 2a),  (A. 2b),  and  (A. 2d),  subject  to  a  given  P(r,  0). 

For  j  *  2,  3,  ...,  N  +  1,  the  solution  on  line  j  is  found  using  solu¬ 
tion  data  on  line  j-1,  the  solution  for  j  =  1  being  used  to  start 
the  recursion.  For  each  j  a  2,  the  relevant  forms  of  GJ  and  H.  (as 
components  of  F),  2  £  i  s  M,  are  therefore  exactly  those  in  Eqs.  (A. 32) 
[cf.  Eqs.  (A. 24)],  variables  with  j-1  suffices  being  known.  From 
Eqs.  (A. 2b),  (A. 2d),  and  (A. 14),  we  also  have 
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u,  .  a  0  ,  T,  .  =  T0  . 

l.J  l.J  2,j 


(A. 34) 


and 


K 


M+l , j - 1 


*  2  (yj  *  yj-l))  ' 


(A. 35) 


Therefore,  for  j  i  2, 


M+l  2 


M1  ,  ,  r-2  ( TM+1 ,  j  +  TM-H,  i-1  ~  2\ 

Ar  (uM+l,j  ~  ^.j  UM+1, j-1  "  UM,j-l)  j  \  2  / 


P  + 

M+l, j-1  + 


f  (t~± —  +  ? — “ — ^  Ay  -  f  (y,  +  y,  ,) 

\  M+l,j  1 M+l, j-1 /  J  3 


Ki 


Vi 


(Tv 


T..  .  +  T„ 


-  x  ) 


2  Ar  M+l, j  M, j  M+l, j-1  Tl.j-l' 
T.  -  2 


-  E 


^!m±lj_1 


M+l, j-1 


)• 


(A. 36) 
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ym,o  "  0  ’ 


'M,l 


-i  + 


2  T  +  T 

M+l,j  M+l , j - 1 


M,0  2  Ar  ’ 


E. 


X-1  2  Ar  TM+l,j  +  T(H-1,J-1 


(A. 40) 


The  detailed  specification  of  G,  H,  and  J  (*  3F/3x)  is  now  com¬ 
plete  for  all  cases.  The  equation  J6  *  -F  is  solved  for  6  using 
Gaussian  elimination  [Isaacson  and  Keller,  1966].  For  efficiency, 
standard  sparse-matrix  methods  are  employed:  only  nonzero  elements 
of  J  [see  Eq.  (A. 19)]  are  stored,  and  operations  involving  zero 
elements  are  automatically  omitted. 


PROGRAM  LISTING 

FORTRAN  listings  of  all  parts  of  the  computer  code  are  now 
presented.  Table  A. 2  lists  the  main  program,  and  Tables  A. 3,  A. A, 
and  A. 5  list  the  subroutines  BCFUNC,  SPRINT,  and  MSWEEP,  respectively. 


LIST  INC.  OF  MAIN  I'ROORAM 


nftfi  ?;>* 


rable  A. 2 — Continued 


I 
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LISTING  OF  SUBROUTINE  BCKUNC 


LISTING  OF  SUBROUTINE  SPRINT 


51  fl'inbrnv  ■'!<!  i»  i  ii'pp, 

52  nnn'jfdt  inn  t  iw'i»  i  (//// 

51  nnnt,?ui  MiiMIt'ul" 

5A  nonhtCI  it’n  F  I  iff,  * 


LISTING  OK  SUBROUTINE  MSWEEP 


ftji  1  f>* P  I 


Table  A. 5 — Continued 


I  «k  f  lOt) 


Table  A. 5 — Continued 


Table  A. 5 — Continued 


NO  lHHOHSiF  ?n  (•‘f'fl.np  StiM><ntif  | •«(  "s«|  ( l' 

STATCmi  »i|  t)i<FFfu:  ,*i>  <•¥*».,  «.  i  <  | 

SINfLF  Phf  I.IS|II»'  F  |  l>  A  I  |  *'f.  *•  I  S,|M|  U|itu|,.i|.  t 


REFERENCES 


Brode,  H.  L.,  D.  A.  Larson,  and  R.  D.  Small,  Time- Dependent  Model  of 
Flows  Generated  by  Large  Area  Fires ,  Pacific-Sierra  Research  Cor¬ 
poration,  Note  483,  July  1982. 

Countryman,  C.  M. ,  PROJECT  FLAMBEAU  ...  An  Investigation  of  Mass 
Fire  '  1964-196? ) ,  Vol.  1,  U.S.  Department  of  Agriculture,  U.S. 

Forest  Service,  Berkeley,  California,  1969. 

DCPA  Attack  Environmental  Manual ,  Chap.  3,  "What  the  Planner  Needs  To 
Know  about  Fire  Ignition  and  Spread,"  U.S.  Defense  Civil  Prepared¬ 
ness  Agency  and  U.S.  Department  of  Defense,  Washington,  D.C.,  June 
1973. 

Isaacson,  E.,  and  H.  B.  Keller,  Analysis  of  Numerical  Methods,  Wiley, 
New  York,  1966. 

Johnson,  L.  E. ,  and  D.  A.  Larson,  Neglected  Effects  in  Nuclear  War¬ 
fare  Simulation,  Pacific-Sierra  Research  Corporation,  Note  438, 

May  1982. 

Keller,  H.  B.,  Numerical  Methods  for  Two-Point  Boundary-Value  Problems, 

Blaisdell,  Waltham,  Massachusetts,  1968. 

Morton,  B.  R.,  G.  I.  Taylor,  and  J.  S.  Turner,  "Turbulent  Gravitational 
Convection  from  Maintained  and  Instantaneous  Sources,"  Proc.  Roy. 

Soc.  A.,  Vol.  24,  1956,  pp.  1-23. 

Murgai,  M.  P.,  and  H.  W.  Emmons,  "Natural  Convection  above  Fires," 

J.  Fluid  Me'  Vol.  8,  1960,  pp.  611-624. 

Nielsen,  H.  J.,  and  L.  N.  Tao,  "The  Fire  Plume  above  a  Large  Free- 
Burning  Fire,"  Tenth  Symposium  (International)  on  Combustion,  the 
Combustion  Institute,  Pittsburgh,  1965,  pp.  965-972. 

Palmer,  T.  Y. ,  "Large  Fire  Winds,  Gases  and  Smoke,"  Atmos.  Environment , 
Vol.  15,  No.  10/11,  1981,  pp.  2079-2090. 

Small,  R.  D.,  D.  A.  Larson,  and  H.  L.  Brode,  Analysis  of  Large  Urban 
Fires,  Pacific-Sierra  Research  Corporation,  Report  1122,  September 
1981. 

Smith,  R.  K. ,  B.  R.  Morton,  and  L.  M.  Leslie,  "The  Role  of  Dynamic 
Pressure  in  Generating  Fire  Wind,"  J.  Fluid  Mech. ,  Vol.  68,  1975, 
pp.  1-19. 


I 


-120- 

Thomas,  P.  h.  "The  cj 

natzonal  Symposium  on  Combustion  !?at“ral  Fires' "  Uinth.  Inter- 

Pittsburgh,  1963,  pp.  844-851'  Combustion  Institute, 

Yokoi,  s..  Study  on  the  P-t'event”  • 

current,  Japanese  Ministry  of  Cons  truer  %  Set  Upuord 

cute.  Report  34,  November^bO.  ’  Bullding  Research  Insti- 


T 


DISTRIBUTION  LIST 


National  Addresses 

Air  Force  Weapons  Laboratory 
Attn:  Civil  Engineering  Division 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117  (1) 

Air  Force  Weapons  Laboratory 
Attn:  SUL  Technical  Library 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117  (1) 

Mr.  Raymond  Alger 
SRI  International 

Menlo  Park,  CA  94025  (1) 

Mr.  Norman  J.  Alvares 
Lawrence  Livermore  National 
Laboratory 

University  of  California 

P.0.  Box  808,  L-442 

Livermore,  CA  94550  (1) 

Assistant  Director 
Energy  and  Natural  Resources 
Office  of  Science  and  Technology 
Policy 

Executive  Office  Building 
Washington,  D.C.  20500  (1) 

Assistant  Secretary  of  the  Army 
(RD&A) 

Attn:  Deputy  ASA  for  RD&S 
The  Pentagon 

Washington,  D.C.  20301  (1) 

Dr.  Jana  Backovskv 
SRI  International 

Menlo  Park,  CA  94025  (1) 

Mr.  A.  P.  Brackebush 
Forest  Fire  Research 
Northern  Forest  Fire  Laboratory 
Missoula,  MT  59801  (1) 

Mr.  Clay  P.  Butler 
SRI  International 

Menlo  Park,  CA  94025  (1) 


Center  for  Planning  and  Research,  Inc. 
Attn:  Document  Library 
2483  East  Bayshore  Road,  Suite  104 
Palo  Alto,  CA  94303  (1) 

Dr.  Craig  Chandler,  Director 
Forest  Fire  4  Atmospheric  Science 
Research 

U.S.  Forest  Service 
Department  of  Agriculture 


Washington,  D.C.  20013  (1) 

Dr.  Conrad  Chester 

Oak  Ridge  National  Laboratory 

P.0.  Box  X 

Oak  Ridge,  TN  37830  (1) 

Chief  of  Engineers 
Department  of  the  Army 
Attn:  DAEN-RDZ-A 

Washington,  D.C.  20314  (1) 

Dr.  William  F.  Christian 
Underwriters'  Laboratories,  Inc. 

333  Pfingsten  Road 

Northbrook,  IL  60062  (1) 

Civil  Engineering  Center/AF/PRECET 
Wright-Patterson  Air  Force  Base, 

OH  45433  (1) 


Dr.  John  Cockayne 
Senior  Scientist 
Science  Applications,  Inc. 

1710  Goodridge  Drive 

McLean,  VA  22102  (1) 

Defense  Technical  Information  Center 
(DTIC) 

Cameron  Station 

Alexandria,  VA  22314  (12) 

Department  of  Defense 

Command  and  Control  Technical  Center 

Attn:  Office  of  the  Director 

The  Pentagon 

Washington,  D.C.  20301  (1) 


Department  of  Energy 
Attn:  Director,  Department  of 
Military  Application 
Washington,  D.C.  20545  (1) 

The  Dikewood  Corporation 
Attn:  Document  Library 
1613  University  Boulevard,  N.E. 
Albuquerque,  NM  87102  (1) 

Director 

Defense  Nuclear  Agency 
Attn:  Michael  Frankel 
Washington,  D.C.  20305  (1) 


Mr.  Marvin  Drake 
Science  Applications,  Inc. 

1250  Prospect  Street 

La  Jolla,  CA  92037  (1) 

Mr.  Donald  Drzewiecki 
Calspan  Corporation 
P.0.  Box  400 

Buffalo,  NY  14225  (1) 

Factory  Mutual  Research  Corporation 
Attn:  Dr.  Ray  Friedman 
1151  Boston-Providence  Turnpike 
Norwood,  MA  02062  (1) 


Director 

Defense  Nuclear  Agency 
Actn:  LtCol.  David  H.  Thomas, 

USAF 

Washington,  D.C.  20305  (1) 

Director 

Lovelace  Foundation 

5200  Gibson  Boulevard,  S.E. 

Albuquerque,  NM  87108  (1) 

Director,  U.S.  Army  Ballistic 
Research  Laboratory 
Attn:  Document  Library 

Aberdeen  Proving  Grounds,  MD 

21005  (1) 

Director,  U.S.  Army  Ballistic 
Research  Laboratory 
Actn:  Mr.  William  Taylor 
Aberdeen  Proving  Grounds,  MD 

21005  (2) 

Director,  U.S.  Army  Engineer 
Waterways  Experiment  Station 
Actn:  Document  Library 
P.0.  Box  631 

Vicksburg,  MS  39180  (1) 

Director,  U.S.  Army  Engineer 

Waterways  Experiment  Station 
Attn:  Mr.  W.  L.  Huff 
P.0.  Box  631 

Vicksburg,  MS  39180  (1) 

Director,  U.S.  Army  Materials  and 

Mechanics  Research  Center 
Attn:  Technical  Library 
Watertown,  MA  02172  (1) 


Federal  Emergency  Management  Agency 
Attn:  Assistant  Associate  Director 
for  Research 

National  Preparedness  Programs 
Directorate 

Washington,  D.C.  20472  (3) 

Federal  Emergency  Management  Agency 
Attn:  David  W.  Bensen 
Office  of  Research 
National  Preparedness  Programs 
Directorate 

Washington,  D.C.  20472  (45) 

Dr.  Francis  E.  Fendell 

Rl/1038 

TRW 

One  Space  Park 

Redondo  Beach,  CA  90278  (1) 

Fire  Research  Library 
National  Bureau  of  Standards 
Building  224,  Room  A-246 
Washington,  D.C.  20234  (1) 

Mr.  Dick  Foster 
SRI  International 
1611  Kent  Street 

Arlington,  VA  22209  (1) 

Dr.  Robert  Fristrom 
Johns  Hopkins  Applied  Physics 
Laboratory 
Johns  Hopkins  Road 

Laurel,  MD  20810  (1) 

Dr.  Matthew  G.  Gibbons 
5424  Lawton  Avenue 

Oakland,  CA  94618  (1) 


Mr.  Edward  L.  Hill 
*  Research  Triangle  Institute 

P.O.  Box  12194 

Research  Triangle  Park, 

NC  27709  (1) 

Dr.  Dennis  Holliday 
R&D  Associates 
P.O.  Box  9695 

Marina  del  Rey,  CA  90291  (1) 

Hudson  Institute 
Attn:  Library 
Quaker  Ridge  Road 

Croton-on-Hudson,  NY  10520  (1) 

Mr.  Pecer  S.  Hughes 

Los  Alamos  Technical  Associates, 

Inc. 

P.O.  Box  410 

Los  Alamos,  NM  87544  (2) 

Mr.  Robert  G.  Hutman 
Nuclear  Test  Engineering 
Division 

Lawrence  Livermore  National 
Laboratory 

University  of  California 
P.O.  Box  808 

Livermore,  CA  94550  (1) 

Professor  A.  Murty  Kanury 
Department  of  Aerospace  and 
Mechanical  Engineering 
University  of  Notre  Dame 
Notre  Dame,  IN  46556  (1) 

Mr.  Kenneth  Kaplan 

#30  White  Plains  Court 

San  Mateo,  CA  94402  (1) 

Mr.  Samuel  Kramer 
National  Bureau  of  Standards 
Building  225,  Room  B-124 
Washington,  D.C.  20234  (1) 

Mr.  Richard  Laurino 
Center  for  Planning  and 
Research,  Inc. 

2483  East  Bavshore  Road, 

Suite  10^ 

Palo  Alto,  CA  94303  (1) 


Dr.  Anatole  Longinow 
I IT  Research  Institute 
10  West  35th  Street 

Chicago,  IL  60616  (1) 

Los  Alamos  Scientific  Laboratory 

Attn:  Document  Library 

Los  Alamos,  NM  87544  (1) 

Mr.  Stanley  B.  Martin 
Stan  Martin  and  Associates 
860  Vista  Drive 

Redwood  City,  CA  94062  (1) 

Dr.  Clarence  R.  Mehl 
Sandia  National  Laboratories 
Division  7112 
P.O.  Box  5800 

Albuquerque,  NM  87185  (1) 

Mr.  Joseph  E.  Minor 

Texas  Technological  College 

Lubbock,  TX  79408  (1) 

Mr.  H.  L.  Murphy 
P.O.  Box  1727 

San  Mateo,  CA  94401  (1) 

National  Council  on  Radiation 
Protection  and  Measurements 
Attn:  Library 
7910  Woodmont  Avenue 

Bethesda,  MD  20014  (1) 

National  Fire  Protection  Association 
Library 

Batterymarch  Park 

Quincy",  MA  02269  (1) 

Oak  Ridge  National  Laboratory 
Attn:  Librarian 
P.O.  Box  X 

Oak  Ridge,  TN  37830  (1) 

Oak  Ridge  National  Laboratory 
Attn:  Emergency  Technology  Division 
Librarian 
P.O.  Box  X 

Oak  Ridge,  TN  37830  (1) 

Dr.  Fred  Offensend 
SRI  International 

Menlo  Park,  CA  94025  (1) 


Mr.  William  Parker 
National  Bureau  of  Standards 
Building  224,  Room  A- 345 
Washington,  D.C.  20234  (1) 

Professor  R.  K.  Pefley 

University  of  Santa  Clara 

Santa  Clara,  CA  95053  (1) 

Mr.  Laurence  Pietrzak 
Mission  Research  Corporation 
P.0.  Drawer  719 

Santa  Barbara,  CA  93102  (1) 

President 

International  Association  of 
Fire  Fighters 
1750  New  York  Avenue,  N.W. 

(3rd  Floor) 

Washington,  D.C.  20006  (1) 

Chief  Robert  G.  Purington 
Lawrence  Livermore  National 
Laboratory 

University  of  California 
P.0.  Box  808,  L-519 

Livermore,  CA  94550  (1) 

The  Rand  Corporation 
Attn:  Document  Library 
1700  Main  Street 

Santa  Monica,  CA  90406  (1) 

Mr.  John  Rempel 

Center  for  Planning  and  Research 
2483  East  Bayshore  Road, 

Suite  104 

Palo  Alto,  CA  94303  (1) 

Dr.  John  Rockett 

National  Bureau  of  Standards 

Center  for  Fire  Research 

Building  224,  Room  B-260 

Washington,  D.C.  20234  (1) 

Mr.  Harvey  r*.  Rvland 
Ryland  Research,  Inc. 

5266  Hollister  Avenue, 

Suite  324 

Santa  Barbara,  CA  93111 


Dr.  Donald  Sachs 
Kaman  Sciences  Corporation 
2001  Jefferson  Davis  Highway 
Arlington,  VA  22202  (1) 

Mr.  Fred  Sauer 

Physics  International  Company 
2700  Merced  Street 

San  Leandro,  CA  94577  (1) 

Dr.  Don  Scheuch 
430  Golden  Oak  Drive 

Portola  Valley,  CA  94025  (1) 

Mr.  Leo  A.  Schmidt 

Institute  for  Defense  Analyses 

Program  Analysis  Division 

1801  N.  Beauregard  Street 

Alexandria,  VA  22311  (1) 

Mrs.  Ruth  W.  Schnider 
Center  for  Planning  and  Research,  Inc. 
2483  East  Bayshore  Road,  Suite  104 
Palo  Alto,  CA  94303  (1) 

Southwest  Research  Institute 
Department  of  Fire  Technology 
P.0.  Drawer  2851D 

San  Antonio,  TX  78284  (1) 

Dr.  Lewis  V.  Spencer 

National  Bureau  of  Standards 

Center  for  Radiation  Research 

Building  245,  Room  C-313 

Washington,  D.C.  20234  (1) 

Mr.  Walmer  E.  S trope 

Center  for  Planning  and  Research, 

Inc. 

5600  Columbia  Pike,  Suite  101 
Bailey's  Crossroads,  VA  22041  (1) 

Technology  &  Management  Consultants 
330  Washington  Street 
Suite  613 

Marina  del  Rey,  CA  90291  (1) 

U.S.  Army  Combined  Arms  Combat 
Development  Activity 
Fort  Leavenworth,  KA  66027 


(1) 


(1) 


U.S.  Army  Training  and  Doctrine 
Command 
Fort  Monroe 

Hampton,  VA  23651  (1) 

U.S.  Forest  Service 
Attn:  Dr.  A.  Broido 
P.0.  Box  245 

Berkeley,  CA  94710  (1) 

Mr .  Thomas  Watermann 

IIT  Research  Institute 

10  West  35th  Street 

Chicago,  IL  60616  (2) 

Mr.  Carl  Wiehle 

Defense  Intelligence  Agency 

Attn:  WDB-4C2 

Washington,  D.C.  20301  (1) 

Dr.  Forman  Williams 
Department  of  the  Aerospace 
and  Engineering  Sciences 
University  of  California  at 
San  Diego 

La  Jolla,  CA  92037  (1) 

Mr.  C.  Wilton 
Scientific  Service,  Inc. 

517  East  Bayshore  Drive 

Redwood  City,  CA  94060  (2) 


International  Addresses 

Almannavamir  Rikisins 

Reykjavik,  Iceland  (1) 

Bundesministerium  des  Innern 
Graurheindorfer  Strasse  198 
5300  Bonn  1 

West  Germany  (1) 

Canadian  Defence  Research  Staff 
Attn:  Dr.  K.  N.  Ackles 
2450  Massachusetts  Ave.,  N.W. 
Washington,  D.C.  20008  (4) 

Civil  Defense  Administration 

Ministry  of  Interior 

Ankara,  Turkey  (1) 


Civil  Emergency  Planning  Directorate 
North  Atlantic  Treaty  Organization 
1110  NATO,  Belgium  (1) 

Directeur  Organisatie 
Bescherming  Bevoling 
Ministry  of  Interior 
Schedeldoekshaven  200 
Postbus  20011 

2500  The  Hague,  Netherlands  (1) 

Directeur  de  la  Protection  Civile 
Minis tere  de  l’Interieur 
36  Rue  J.  3.  Esch 

Grand-Duche  de  Luxembourg  (1) 

Direction  de  la  Securite  Civile 
Ministere  de  l'Interieur 
18  Rue  Ernest  Cognac 

92  Levallois  (Paris)  France  (1) 

Director 

Civilforsvarsstyrelsen 
Stockholmsgade  27 
2100  Copenhagen  0 

Denmark  (1) 

Forest  Fire  Research  Institute 

The  Information  Center 

331  Cooper  Street 

Ottawa,  Ontario  KIA  043 

Canada  (1) 

Brigadier  I.G.C.  Gilmore 
Director,  Australian  Counter 
Disaster  College 
Mount  Macedon,  Victoria  3441 
Australia  (1) 

The  Head  of  Sivilforsvaret 
Sandakerveien  12 
Postboks  8136 
Oslo  dep 

Oslo  1,  Norway  (1) 

Home  Office 

Scientific  Advisory  Branch 
Horseferry  House 
Dean  Ryle  Street 
London  SW1P  2AW 

England  (1) 


Jefe,  Section  de  Estudios  y 
Planification 
c/Evaristo  San  Miguel,  8 
Madrid-8 

Spain  (1) 

Mlniscero  dell  Incerno 
Direzione  Generale  della 
Protezione  Civile 

00100  Rome,  Italy  (1) 

Ministry  of  Social  Services 
11  Spartis  Street 

Athens,  Greece  (1) 

Secretaire  d' Administration 
Ministere  de  l'Interieur 
Direction  Generale  de  la 
Protection  Civile 
Rue  de  Louvain,  1 

1000  Brussels,  Belgium  (1) 

Dr.  Ing.  P.G.  Seeger 
Forschungsstelle  fur 
Brands  chut z  t echnik 
University  of  Karlsruhe  (TH) 

75  Karlsruhe  21 
Postfach  63380 

West  Germany  (1) 

Servico  National  de 
Proteccao  Civil 
Rua  Bela  Vista  a  Lapa,  37 
1200  Lisbon,  Portugal  (1) 

Dr.  Vilhelm  Sjolin 
Director  of  BRANDFORSK 
The  Swedish  Fire  Research  Board 
S-115  87  Stockholm 

Sweden  (1) 

Stato  Maggiore  Difesa  Civile 
Centro  Studi  Difesa  Civile 
Rome,  Italy 


(1) 


PSR  Report  1210 


ANALYSIS  OF  THE  LARGE  URBAN  FIRE  ENVIRONMENT 
Part  II.  Parametric  Analysis  and  Model  City  Simulations 
Summary 


By 

D.  A.  Larson 
R.  D.  Small 


November  1982 


Final  Report 

Contract  EMW-C-0747,  Work  Unit  2564E 


For 

Federal  Emergency  Management  Agency 
National  Preparedness  Programs 
Washington,  D.C.  20472 


FEMA  Review  Notice 

This  report  has  been  reviewed  in  the  Federal  Emergency  Management  Agency 
and  approved  for  publication.  Approval  does  not  signify  that  the  contents 
necessarily  reflect  the  views  and  policies  of  the  Federal  Emergency  Management 
Agency. 

Approved  for  Public  Release:  Distribution  Unlimited 


PACIFIC-SIERRA  RESEARCH  CORP. 

12340  Santa  Monica  Blvd.  •  Lot  Angeles,  CA  90025  •  (213)  820-2200 


S-l 


SUMMARY 


This  report  considers  the  large-fire  environment  that  would  occur 
in  an  urban  area  subject  to  a  nuclear  weapon  explosion.  The  effects 
of  system  parameters  are  explored  in  a  sensitivity  study,  and  results 
for  three  model  cities  are  presented. 

The  examples  are  characterized  by  extensive  areas  simultaneously 
burning,  strong  buoyancy,  and  large  temperature  gradients.  Several 
such  fires  occurred  during  World  War  II.  Though  those  fires  were 
dramatic  in  intensity  and  destructiveness,  each  involved  a  relatively 
small  area.  A  nuclear  weapon  explosion  could  generate  a  far  larger 
area  fire  and  a  more  severe  fire  environment.  This  report  is  intended 
to  define  such  large  fires.  The  results  should  be  applicable  for 
damage  evaluation,  formulation  of  shelter  requirements,  and  rescue 
planning. 

The  calculations  are  based  on  the  theory  developed  in  Part  1  of 
this  report,  which  is  applicable  to  the  fire  zone  and  the  volume 
immediately  above  it  (turning  region).  The  effects  of  variable  area 
heating,  turbulence,  strong  buoyancy,  large  temperature  changes,  and 
radiation  are  treated.  The  induced  fire  winds  and  rapid  temperature 
changes  at  the  fire  periphery  are  uniquely  determined  by  the  use  of 
jump  conditions.  Simulations  of  the  Hamburg  firestorm  and  a  large 
Flambeau  fire  agreed  well  with  available  data. 

The  parametric  analysis  considers  a  large  area  fire  and  the  ef¬ 
fects  of  fire  size,  heating  rates,  mixing  coefficients,  and  hot  gas/ 
smoke  radiation.  The  results  show  the  influence  of  those  variables 
on  the  induced  fire  winds,  mean  temperature,  and  pressure  gradients. 

In  general,  an  increase  in  either  the  fire  size  or  heating  rate  raises 
the  mean  temperature  levels  and  the  induced  fire-wind  velocities. 

For  the  larger  heat  release  rates  or  fire  sizes,  the  attendant  in¬ 
creases  in  mean  temperature  and  velocity  are  limited  by  compressi¬ 
bility  effects. 

Fires  such  as  may  result  from  a  megaton-yield  explosion  are 
analyzed  for  three  model  urban  areas.  Each  city  is  characterized  by 


a  high-density  center,  a  surrounding  belt  of  mixed  residential/ 
industrial  construction,  and  a  lower  density  suburban  belt.  Each 
model  city  portrays  a  different  degree  of  building  density  and  urban 
sprawl.  The  results  illustrate  how  a  particular  city  geometry  affects 
the  velocity  and  temperature  fields  for  a  given  fire.  An  additional 
series  of  computations  considers  reduction  of  the  fire  area  by  severe 
blast  damage  and  debris  formation.  For  those  calculations,  complete 
burning  was  allowed  in  an  annular  area,  with  the  fire  intensity  signif¬ 
icantly  reduced  in  the  center. 

Finally,  the  model  is  employed  to  estimate  the  behavior  of  the 
velocity  and  temperature  fields  as  a  function  of  fire  evolution. 

Those  calculations  may  indicate  the  most  appropriate  periods  for 
effecting  rescue  operations  as  well  as  provide  an  estimate  of  time- 
dependent  shelter  loadings. 


